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Expressions for decay amplitudes of B°,B~ and B® mesons to two light pseudoscalar or vector 
mesons, including isosinglet mesons r),r)',u),(f>, are obtained using Soft Collinear Effective Theory 
at LO in l/rat- These are then used to predict unmeasured branching ratios, direct and indirect 
CP asymmetries in B°,B~ and B® decays into two light pseudoscalars, following a determination 
of nonperturbative SCET parameters from existing data using a y^-Ht. A separate discussion of 
indirect CP asymmetries in penguin dominated B° — » j/ ' Ks,l,^° Ks,l decays is given. 



I. INTRODUCTION 

Hadronic two-body B decays constitute both a probe 
of the electroweak structure of the Standard Model as 
well as a testing ground for our understanding of QCD 
dynamics. In B — ► MM' decays with two energetic 
light final mesons M, M' going back-to-back in oppo- 
site directions, the amplitude factorizes in the limit 
TO;, 3> Aqcd, ttim,m'- This had made possible the great 
advancements in our theoretical understanding of B — > 
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The observation of factorization initially relied on a 2- 
loop proof (for earlier attempts see @). Following the 
advent of Soft Collinear Effective Theory (SCET) jf the 
factorization was then shown to hold to all orders in as 
at leading order in 1 /nib 0, IE IE IE , with a possible 
exception of the contributions from intermediate charm 
quark states 0, EE [H EH . 

In this paper we present the first SCET analysis of two 
body B decays where one of the final mesons contains 
an isospin singlet admixture to the wave function. This 
includes all the decays with rj, 7/ or <j), uj mesons in the 
final state. As we will show the factorization in the limit 
To;, — > 00 which was obtained for nonisosinglet final states 

111 Refs. 01 una, generalizes to the case of isospin 

singlet final states. The expression for the B — > MM' 
amplitude is then schematically 
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with £g) denoting the convolutions over momenta frac- 
tions, 4>m the light-cone distribution amplitude (LCDA) 
of a meson that does not absorb the spectator quark, 
C, BM the soft overlap function, Cj M the function de- 
scribing the completely factorizable contribution, and 
T^.j the corresponding hard scattering kernels, while the 
ellipsis denotes the charming penguin contributions. 

Eq. already contains additional contributions aris- 
ing from purely gluonic configurations that are specific to 



decays into pseudoscalar isosinglet states. These for in- 
stance lead to new jet functions in the SCETi — >SCETn 
matching in addition to the ones found in 0, so that 
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both C Bv ' and C,j'"' ' 

receive leading order gluonic con- 
tributions. The phenomenological importance of rj glu- 
onic content was emphasized by a number of authors 
[24l I2H l2t| after a discovery of surprisingly large branch- 
ing ratio ~ 10 -4 for semi-inclusive B — > vj ' X decays some 
years ago [23]. These expectations are at least partially 
confirmed by our SCET analysis. The gluonic contribu- 
tions to the amplitude where the spectator light quark 
in the B meson is annihilated in the weak vertex are of 
leading order in the l/m& and as(rrib) expansions. This 
result also confirms the discussion of gluonic contribu- 
tions presented in the QCD factorization calculation of 
B — > rj' K modes |0, Q|, where a similar contribution 
was taken into account as part of B —> rj' form factor. 
Unlike the authors of Ref . [TH . fl4| , which were forced 
to assign a rather arbitrary size of either 0% or ~ 40% 
for the gluonic contribution to the B — > -q' form fac- 
tor, we are able to use the wealth of new data from the 
B-factories and fit for the corresponding SCET nonper- 
turbative parameters. 

Following we do not expand the jet functions in 
terms of as(vAmb) and treat the corresponding func- 
tions as nonperturbative parameters that are determined 
from data al ong with the charming penguins. Although 
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1 1 5| have argued that the charming penguins 
are perturbatively calculable, the more conservative ap- 
proach of determining the charming penguins from data 
avoids any inherent uncertainties associated with possi- 
bly leading order long distance effects [EElElliS- Tak- 
ing a conservative approach is especially important, if one 
aims at interpreting hints of beyond the Standard Model 
contributions to the observables in AS — 1 processes: the 
nK puzzle [2EI2EI2E] an d the deviations of S parameters 
in penguin dominated modes Ks rj ', K ^tt , Ks<t> from the 
naive expectation of S ~ sin 2/3 |3l|, |33, y3| . We devote 
the second part of our paper to these phenomenological 
considerations. In the numerical studies we assume that 



2 



1/nib corrections are of typical size ~ 20% 

In this paper we provide expressions for B — > PP, PV 
and decays at LO in the 1/rrib and as(jrib) expan- 
sions, while in the numerical estimates we restrict the 
analysis to the decays into two pseudoscalars. Similar 
numerical analyses of B — > tttt, ttK and B — > if if decays 
using SCET were performed previously in Refs. |9l Il0|. 
In addition to these modes we also discuss decays into 
final states with rj and rf mesons, usin g a X 2 -fit to ex- 
tract the SCET parameters. Since in a subset of 
measured observables rather than a x 2_n t was used to fix 
the values of SCET parameters, the numerical results of 
the two analyses do not match exactly, but they do agree 
within the errors quoted. Furthermore SU(3) flavor sym- 
metry for the SCET parameters is needed at present due 
to limited data for decays into 77 and r\' mesons. 

Assuming only isospin symmetry, the amplitudes for 
AS* = decays with isosinglets in the final states, i.e. 
the decays B —> irrf ' and B — > rf>rf\ are written 
in terms of eight new real nonperturbative SCET pa- 
rameters beyond those describing the B — > tttt system 
(see Eqs. Q43 [l -Q49 [l below) at leading order in 1/mb- In 
total there are nineteen observables in B — > nrf ' and 
B —* rf'rf} decays, only four of which have been mea- 
sured so far. Therefore, we are forced to use SU(3) flavor 
symmetry to reduce the number of unknowns. The situ- 
ation is very similar for AS* = lS-> Krp •* decays with 
a total of ten observables, seven of which have been mea- 
sured so far, while the amplitudes are expressed in terms 
of eight new real parameters beyond the ones already 
present in B — > tttt and B — > ttK decays. Using SU(3) 
flavor symmetry greatly reduces the number of parame- 
ters. In this limit all the B — > PP decays are described 
at LO in l/m& in terms of only eight real unknowns 
(four without isosinglet final states). This should be con- 
trasted with the conventional use of SU(3) decomposition 
that leads to eighteen reduced matrix elements (nine for 
each of the two independent CKM elements combination, 
see Appendix [UJ and therefore to 35 real unknowns and 
one unobservable overall phase. We thus choose to work 
in the SU(3) flavor limit and perform an analysis of all 
presently available data on B — ► PP decays from which 
we predict the values of yet unmeasured observables in 
both AS = and AS = 1 decays of B°, B~ and B° s . 

The paper is organized as follows: in Section [H] the 
matching of QCD^ SCETi — ►SCETn and the resulting 
amplitudes for two-body decays into light pseudoscalar 
and vector mesons (including decays into flavor singlet 
mesons) at leading order in 1/mt are derived. The phe- 
nomenological implications of these results for B~,B° 
and P|? decays into two light pseudoscalar mesons are 
then developed in Section ILTTl Notation used throughout 
the paper is collected in Appendix lAl while Appendix 151 
deals with the Dirac structure of the operators multiply- 
ing jet functions. Finally, in Appendix ICl our results are 
rewritten in terms of SU(3) reduced matrix elements and 
a translation to the diagrammatic notation is made. 



II. TWO-BODY B DECAYS IN SCET 

The starting point is the effective weak Hamiltonian at 
the scale fi ~ m& for the AS = 1 two body B decays [34[ 

io,7 7 ,8 9 

v p— u,c i—Z 
(s) 

where the CKM factors are X p — V p bV* and the stan- 
dard basis of four-quark operators is 

°i = (P b )(sp)-, Ol = (ppb a )(s a pp)-, 
03,5 = (sb)(qq) T , 4 ,6 = {s a bf)){qpq a ) T , ^ 

07,9 = ^Y(sb)(qq)±, 08,10 = ^(SaM^aU, 

with the abbreviation (gi7 AI (l-75)g2)(g37 M (lT75)94) = 
(9i ( ?2)(93<?4)=f- The color indices a,f3 are displayed only 
when the sum is over fields in different brackets. In the 
definition of the penguin operators 03- 10 in Q there is 
also an implicit sum over q — {u, d, s, c, b}. The electro- 
magnetic and chromomagnetic operators are 

°{77,8 S } - -^Sa^ieF^, gG^}P R b. (4) 

The weak Hamiltonian for AS — decays is obtained 
from 10-0 through the replacement s — > d. We will 
be working to leading order in asijrib) so the values for 
the Wilson coefficients in (|2J| are given at LL order in 
the NDR scheme for a s {m z ) = 0.119, a om = 1/128, 
m t — 174.3, even though NLL values are available [34| . 
At the scale /a = nib = 4.8 GeV the Wilson coefficients 
Ci for tree and QCD penguin operators are j33| 

Ci- 6 (m b ) ={1.110, -0.253, 0.011, -0.026, 0.008, -0.032}, 

(5) 

while for electroweak penguin (EWP) operators [HH 

C 7 -w(m b ) = {0.09, 0.24, -10.3, 2.2} x 10" 3 , (6) 

and for the magnetic operators 07 7 (?rih) = —0.315, 
C 8g (m b ) = -0.149. 

The effective weak Hamiltonian of full QCD is 
matched to the corresponding weak Hamiltonian in 
SCET. In the two-body B -> M X M 2 decays there are 
three distinct scales, the hard scale ~ mt due to the en- 
ergy available to the decay products in the B rest frame, 
the typical hadronic soft scale A, and the hard collinear 
scale y/rribA. This last scale corresponds to the typi- 
cal momentum transfer needed to boost the spectator 
quark in the B meson with soft momentum k ~ A so 
that it ends up in the final meson M\ with a momentum 
p ~ mf,(A 2 ,l,A), where A = A/m&. The notation here 
is p^ = (n ■ p,n ■ p,p±), with the meson Mi containing 
the spectator quark going in the = (1,0,0,-1) di- 
rection, while Mi goes in the opposite — (1,0,0,1) 
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direction. The presence of three distinct scales leads to 
a two-step matching procedure [35| • First hard scale raj, 
is integrated out so that the effective weak Hamiltonian 
© in full QCD is matched to the Hamiltonian in SCETi 
where the scaling of jet momenta in n and n directions 
are m b (X, and mb(l, \, respectively. This is 

then matched to a set of nonlocal operators in SCETn 
with jet functions as Wilson coefficients |35j. In the fol- 
lowing two subsections we closely follow the matching at 
LO in l/m b performed in Ref. 9J, extending it at the 
same time to a larger set of operators including the oper- 
ators that contribute only for the isosinglet final states. 

A. Matching to SCETi 




FIG. 1: The configuration with intermediate on-shell charm 
quarks annihilating into two collinear quarks going in the op- 
posite directions can lead to nonperturbative long distance 
contributions due to soft gluon exchanges [jj- The double 
line denotes heavy quark, spectator quark is not shown. 



To have a complete set of leading order contributions 
in the second step, we keep in the SCET/ Hamiltonian 



9] 



V2 



n,n i 



(7) 



both the leading order operators Qf) as well as a subset 



'if 

of relevant subleading operators Q^j in the VA expan- 
sion, where f = d (f = s) for AS — (AS* = 1) pro- 
cesses. In (J7J) the charming penguin contributions |23j 
are isolated in the operator Q 5c - Since 2m c ~ m b the in- 
termediate charm quarks annihilating into two collinear 
quarks (see Fig. ^) is a configuration where the interme- 
diate on-shell charm quarks have a small relative veloc- 
ity v and can lead to long distance nonperturbative ef- 
fects due to the exchange of soft gluons. These contribu- 
tions are as{2m c ) f {2m c / mb)v parametrically suppressed 
|il Eil ED < w ith /(2TO c /mb) a factor encoding that only 
in part of the phase space the charm quarks have small 
relative velocities. The view of BBNS [J, [l5|| is that this 
phase space suppression of the threshold region is strong 
enough so that nonperturbative contributions are sub- 
leading. Bauer et al. on the other hand argue 
that since 2m c /m& ~ 0(1) then also f(2m c /mb) ~ 1. 
Sizable charming penguin contributions have also been 
found in recent light-cone calculation [3(j . The most con- 
servative approach is to introduce new unknown param- 
eters describing charming penguin contributions that are 
then fit from data, so this is the approach we will follow. 

In J7J) the same notation as in Refs. jjj 0] is used. 
The set of leading order operators in this notation is 



O {0) - Is 



Qhlfis = [Sn.wi fiPLb v ][qn,u 2 ^PL,R.qn,u 3 }, 



(8) 



and the gluonic operator 

QfJ = m b [s nm flP L b v } Tr[B^ 2 B^ 3 ]ie ±flV) (9) 

where the trace is over color indices. Above, an implicit 
summation over light quark flavors q — {u, d, s} is under- 
stood. The operators Q^ 6Stgs contribute only to isospin 
singlet mesons and were not needed in 0, 0] . The quark 
fields in Equations © and © already contain the Wilson 
line together with the collinear quark field 



q n .u = [S(u - fi ■ V)W$& 



(?)i 



(10) 



with the usual bracket prescription that V operates only 
inside the square brackets |2fl. We also define a purely 
gluonic operator related to the (n, _L) component of the 
gluon field strength 



'—'in f i J 



[W^in ■ D c , n ,iD^ x ]W n S(u> - n ■ V^)} . 



(11) 

In © the operators with T A ® T A color structure were 
not listed, since they do not contribute to color singlet 
final states. Similarly, a gluonic operator with e_\_^ v — > 
5_L/w in EJ is n °t considered, as it does not contribute 
to P, V final states due to parity. The Q^) operators for 
the AS = weak Hamiltonian are obtained from (JSJ) , Q 
with the replacement s — > d. 

The relevant 0{\f\) operators are 



Qu = — W^^gBi P L b v \\s^ 



LUn, 



QL,3 S = — [s^ujjg ^P L b v ][u n 



m b 



iAs - —[Qn,uJgB ntUi PLb v ][Sn,oj 2 '. 



J Lq 



n, uj 3 \, 



Q5s.6s = [Sn^Jg Pn,UiPLbv\[qn,u 2 ¥- P L Mn, 0,3], 

Q7J = — [un^iigB^PLKWsn,^ 
_2 

— [q n , ujgB^PLbv] [s n , U2 i^P R q 



(12) 
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FIG. 2: Matching to Q g J vanishes at LO in l/m b with a 
representative diagram given in a). In b) nonzero diagrams 
that contribute at LO in 1/mj to Q^J matching are given. 
The square represents insertion of 0$ g . The up- (down-) going 
wiggly-solid lines denote n (n) collinear gluons. The offshell 
gluon (wiggly line) in the third diagram can also attach to 
any of the cross vertices. 



with a Wilson coefficient that is a function of 2m c /mb 
and can be found from results in [l3j|. Similarly, the 
remaining diagrams of Fig. Olwith one n gluon emitted 
from b or n collinear quark line are part of the matching 
onto at least V~\ suppressed operators. 

This leaves us with the set of operators and H12|) 
that contribute at LO in 1 /mb and as (nib). Note that 
the Wilson coefficients in (JTJl already contain the CKM 
elements. Furthermore, the basis of operators JSJl and 
(|12|l is a minimal choice, where the operator relation 
O g )10 = §[0^ + O c 21 ] - |0 3) 4, valid after off-shell b 
quarks are integrated out, was utilized. This leads to the 
following tree level matching for the four-quark operators 




FIG. 3: The power suppressed contributions with off-shell 
charm quarks emitting collinear gluons in n direction. 
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(15) 



while NLO matching can be obtained from Q, Q . The 
tree level matching of the dimension seven operators 
leads to 



and the operator with additional n collinear gluon fields 

Q g X ) - ~2[s n ^g^ 4 P L b v ] Tr[B^ 2 B^]ie ±IM/ , 

(13) 

while, quite similarly to ©, we do not consider the ad- 
ditional gluonic operator with g± (ll , instead of e_i_^ v in 
(II 3p , since it does not contribute to P, V final states due 

to parity conservation. As before are obtained by 
making the replacement s — > d in the above definitions. 

The Wilson coefficients cf' for the operators in (JSJ 
and 4s2«^ or the operators in i|12|) are known to NLO in 
as(mi) 0,0123' while the rest of the Wilson coefficients 
b\ for the operators in l|12|l are at present known only to 
leading order. In the phenomenological part of our paper, 
Section IlIII we work at LO in the as{mb) expansion. At 
this order there are no contributions from QgJ and Q g ^J 
operators (JHJ, i|13[) since the tree level matching shown 
in Fig. |21 and given explicitly below, is already as(mb) 
suppressed. 

We also work at leading order in A/mb- In particular 
this means that we systematically neglect the contribu- 
tions to B — ► Prfi > amplitudes with two collinear gluons 
being emitted from the charm loop (see Fig. Ola)). Here 
the intermediate charm quarks are off-shell, unlike the 
charming penguin configuration, Fig. ^ The diagram in 
Fig. |2| a) matches onto the (VX) 2 suppressed operator 
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(16) 



while the NLO contributions from tree operators only 
was recently obtained in |37j . 

For completeness we also list the result of tree level 
matching for Q g °J and QgJ operators ©, following 
from the diagrams in Fig. [21 



[s n , ui fiP L b v ] Tr[S^ 2 m • D n B£"]e^ u , (14) 



-(/) 



0. 



(17) 
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and 



and 



(/)^ as{m b ) 



2 1- 



16C F 
1 



2 + z 



TV 2 / (1 - zu){l - zm) 



(18) 



where z = wi/toj, and u — Lo^/nib, with a>i,3 the fraction 
momenta in (0, while C F = (N 2 - l)/2N and N = 3. 
Even though the matching result (|17|l , l|18|) , which is al- 
ready NLO in as (mi,), will not be needed for our phe- 
nomenological discussions in Sect ion ITTT1 it does represent 
first nonzero contribution from 0% g weak operator to the 
matching. Since this chromomagnetic operator may be 
enhanced in new physics models (in MSSM this is due to 
the fact that the chirality flip can be performed on gluino 
line instead on the b line HE 13), we provide the match- 
ing as an important input for future studies of such new 
physics effects in two body B decays using SCET. 



B. Matching to SCETn 

The matching of SCETi onto SCETn is performed by 
integrating out the degrees of freedom with p 2 ~ Amj,. 
To do so it is useful to perform a redefinition of fields 
£fi — > lfi£n, Afi — ► YftAfiYl, where n denotes either 
n or n direction, while Y n are Wilson lines of n ■ A us 
SCETi usoft gluon fields. After the redefinition the usoft 
fields decouple from collinear fields in the leading order 
SCETi Lagrangian 

HIE El They sti11 

appear in the op- 
erators Qi s G {Qi]i ■ ■ ■ > Q& s \ Qgs j Qis i ■ ■ ■ > Qs s > Qgs}, 
but only in the combination Y%b v . Thus n collinear and n 
collinear sectors decouple [j| and the operators Qis factor 
into 1 [| 



Qis — QisQis- 



(19) 



Since the n and n operators do not communicate through 
usoft gluons, we can focus only on the matching of Q" s op- 
erators onto SCET n operators. The soft spectator quark 
in the B meson is boosted to a collinear quark in the 
n— direction through an application of the 0(\fX) sup- 
pressed Lagrangian [sol ] 



(20) 



in the T products (see also Fig. 0} 



T M/ = / d?y d 4 y' T[Qf\0),iC^ n (y') + i£$(y'), 



2,i/ 



(22) 



The explicit forms of subleading SCET Lagrangians 



C 



(i) 



dig , can be found in [3£ 



Note that all the 



terms in T\^f and ?2,i/ have the same y/X scaling. The 

superficially enhanced T product T[Q^j(0),iC^ (y)] 
gets suppressed once the external momenta are restricted 
to SCT^Tn scaling, since it involves an odd number of 

The matrix elements (M\T\jf\B) lead to endpoint sin- 
gularities, signaling a presence of soft-ove rlap contribu- 
tions both for nonisosinglet final states [35l l4fl ItH l42l IIE 
13 EE HE HI as well as for isosinglet M. Following [!j 
we therefore define the matrix elements of T\^f as new 
nonperturbative functions to be determined from data 



(M\T 1 , if \B) = C l B f M ( BM . 



(23) 



The final meson M can be either P, or Vfi , while 
(V±\Tijf\B) vanish at leading order 0. The coefficients 
Cf^ M = 1, ±l/\/2 describe the flavor content of the final 
meson. 

To facilitate the calculation of T 2 ^ / the Q™j^ operators 
are stripped of the soft quark fields (^f a Y 1 \b v )' la ((Y^[b v )' la 
in the case of Q™^) and similarly C9^ is stripped of 
(q' us Y n y b . The common matching result is then (see ap- 
pendix [5J| 

T[(^ ) W„)^i 5 <LP fl , i ] ia (0)[^^Wt e ( 3 ')]^ (y) = 
iS ab 5(y + )S^(y ± ) f dx J ^e*+»-/ a 
^[^J{z,x,k + ){yiP L ^'l) :ll [q ntXU] flP L ,R 



Lj ± ( Z , x , k+ )^p R , L1 z^y i [q,, 



IB In, — XUJj 



i£ £M + / dSTlQ^iO), iCg*\y)], J g (z, x, k+) and J' g (z, x 



+S q>q> J g (z, x,k+)( iP L ,Rll) 31 ^\Bi%JBi% J 

+5 q , q ,J' g (z, x, k + )(^P L ,Rllf ItPiU ' B i^} + ' ' • ' 

(24) 

where the trace is over color indices, while ellipses denote 
operators with nontrivial color structure that do not con- 
tribute to color singlet initial state. This generalizes the 
result of Ref. [jj to the case of isosinglet final state, 
where also the operators with two interpolating gluon 
fields contribute, leading to two additional jet functions 



(21) 



1 Operators > Qis ,8s factorize into a sum ^2 q Qq,i s Qq,is > but 
we suppress this sum in the notation. 



). At tree level J(z,x,k + ) = 
J±(z,x,k+) = 5(z — x)as^C f /(N r xk+), with 1-loop 
corrections calculated in [42L |46j, and J g (z,x,k + ) = 
S(z - x)a s 2n/(N c k+), J' g (z, x, k+) = 0. 

The matching of different operators is ob- 

tained from by multiplying with (PR^ a Y^b v ) la (by 

(P L Y^b v ) ia in the case of Q^L) and with (q' us Y n) jb - The 







FIG. 4: The T-products of weak Hamiltonian operators (large black dot) with C^„ i!2( )t that boosts soft quark (dashed line) to 
n collinear quark (blue/dark gray solid line). Two representatives of soft overlap contribution diagrams (12 111 with a) outgoing n 
collinear quarks and b) outgoing n collinear gluons are shown. Diagrams c) and d) give factorizable contributions corresponding 
to 1221 1. The red/light gray solid lines denote n collinear quarks, the blue/dark gray wiggly and solid lines the n collinear gluons, 
while double lines denote heavy quarks. The diagrams b) and d) contribute only to r/, rj going in n direction. The as(mb) 
suppressed diagrams with n collinear gluons are not shown. 



final result for B — > M\Mi amplitude depends on only 
two jet functions, J and J g . The operators multiply- 
ing J' g do not contribute since the matrix elements of 
Tr[B„ xul ■ B^xui] between vacuum and pseudoscalar me- 
son or between vacuum and vector meson vanish due to 
parity conservation. That Jj_ jet function does not ap- 
pear in B — > M\Mi amplitudes when are nonisos- 
inglet mesons was already shown in 9], and is true also 
when either of the two outgoing light mesons is an 
isosinglet meson. Namely, the operator multiplying Jj_ 

cancels in the matching of <2ij ■ ■ ■ Qei > Qgs operators 
because 7"7^7j_ q = 0. It also cancels in the matching of 



the QjsSs operators, since ^Pt7"7^ = {g°f 
leads to a vanishing term 



a/3 



(9 



a/3 



a/3 



)[Qn 



^7a«ft] 



0. 



(25) 





FIG. 5: The SCETi — >SCETn matching onto operators con- 
tributing only to B* decays with isosinglet mesons in n di- 
rection. The weak SCETi operator arising at tree level from 
Os g insertion (large black dot) is enhanced due to one less n 
gluon field compared to operators in Fig. |1J which is compen- 
sated by additional soft gluon in the initial state. Additional 
diagrams similar to Fig. 0]b), d) are possible for isosinglet 
final state in n direction. 



where for the last equality the relation 7"7]_7_l q = was 
used once again. Furthermore, matching of the Q™f op- 
erators leads to a soft operator [q U s^PR"ifb v ] multiplying 
J, J g or J'g jet functions. The matrix element of this op- 
erator between vacuum and pseudoscalar \B) states van- 



ishes. Thus Qys8 S operators contribute only to B* 
not at all to B — ► M\M% decays. 

We note on passing that B* decays into two light 
mesons would also receive, at leading order in l/m&, con- 
tributions from non-valence Fock state of the B meson 
with an additional soft gluon, i.e. from the configurations 
with incoming heavy quark, soft spectator quark and soft 
gluon. The SCETi ->SCET n matching is shown on Fig. 
03 Since the weak SCETi operators are 1/VA enhanced 
compared to Ti,i/,T2,j/ (the T products in n direction 
are the same as in (I21II . (|22p. while in n direction the 
weak operator has only one B^ field) , this allows for an- 
other external A us field. Due to Dirac structure these 
operators do not contribute to B decays, quite similar to 
equivalent contributions found in B* — > P7 with P 



and 



an isosinglet. Note that these contributions would not be 
present in B* — » M form factors. 

We are now ready to write down the final result for 
the B -> MiM 2 amplitude at LO in l/m b in SCET that 
will extend the result of Ref |9j to the decays including 
isosinglet mesons. This result exhibits two levels of fac- 
torization: the n direction collinear modes decouple from 
the rest, and for the 2~2,i/ term the n collinear fields de- 
couple from soft fields. Because of the factorization, the 
amplitude has a very simple form Q expressible in terms 
of several universal nonperturbative functions such as the 
£ SM it2*5|) function and the twist-2 light-cone distribution 
amplitudes (LCDA) that are for pseudoscalar final 
state defined through (see e.g. pSOjl 

{P a {p)\q n ^X a M n ^ 2 \Q) = -if P 2E0 pa (u), (26) 

with M = X a P a ' a 3 x 3 matrix of light pseudoscalar fields 
P a given in I|C1|) . and for the longitudinally polarized 
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vector meson final state through 

(V((p)\q n ^X a fiq' n ^ 2 \0) = lf v 2E4> V Hu) 



(27) 



Amplitudes for decays to transversely polarized vector 
mesons do not receive leading order contributions from 
operators (jHJ, lfl2*|> (see also discussion around Eq. J2HJ). 
For the gluonic operator on the other hand 



ie x ^{P(p)\Tr{B^ ul B^_ U2 }\0) = -^fC^fp^u), 

_ _ (28) 

where /* = yfflZf^, /* = f v ,/VS and fj, = for 
other pseudoscalars. Here f Vq and f Vs are the decay con- 
stants corresponding to qq = (uu + dd)/y/2 and ss axial 
currents in l|26|) respectively. The matrix elements of glu- 
onic operators between vacuum and vector meson states 
vanish. In all of the above definitions the integration 
Jg 1 du5(uJi — un ■ p)S(ui2 + uh ■ p) on the r.h.s. is implic- 
itly assumed. The LCDA $p(tt) is related to the LCDA 
(j)p g {u) used in |5l| through 



u(l — u) 



(29) 



and coincides with the definition in [53 ■ It has the sym- 
metry property 4*p( u ) = — 0p(l — u ); as can be easily 
seen from (|28f) . Finally, the relevant B meson LCDA is 
11 



(30) 



Combining the SCETi -> SCET n matching {2H with 
the definitions of the nonperturbative matrix elements 
and using the relation 



(0\(q s Y n )\ x _^(l- l5 )(YX)\o\B) = 

-if B m B [ dr+e- tr + x -/ 2 4+(r + ). 



(31) 



for the gluonic operator in (|24|) . we arrive at the expres- 
sion for the decay amplitude into pseudoscalar or longi- 
tudinally polarized vector mesons valid at leading 
order in 1/rrif, and to all orders in as(mt) 



A(B -> M X M 2 ) = 

~J^ m B {fMi J du(j) Ml (u) J dzTu{u, z)C BM2 (z) 
+ f Ml Jdu0MAu)T lc (u)( BAh 
+ fl h J du¥ Mi (u) J dzTt,{u, z)( BM * (z) 
+ fl h Jdu¥ Mi {u)T 3 K {u)Q BM * 



(32) 



where A^, lM2 is a term denoting the nonperturbative 
charming penguin contributions to be discussed in more 
detail below, see Eqs. l(i7 )l -l|^i | . ( BM was defined in 
while ( BM (z) is given in terms of the jet functions 124|) 
as 



/-BM 1 



fsf. 



M 



m B 



dk+cj)+(k + ) C BM dx4> M (x)J(z, x, k+) 



r<BM x 

C 9 4 



^ I <l.r<I> ! \ l (.r)Jj :..)■./,■_ ! 



(33) 

for both isosinglet and nonisosinglet mesons M. The 
term with gluonic jet function J g (z,x,k + ) contributes 
only when M = r)q,T) s . The coefficients Cf^' 1 parame- 
terizing the relative weights of the two contributions in 

Cf M (z) are 



BM 


C BM 


(jBM 


B%,B-r) q 


1 


2 


B°rj s ,B-r) s 





1 


B% 


1 


1 


B% 





2 



In all other cases, when M is either an isospin nonsinglet 
pseudoscalar meson or a vector meson, we have simply 
C BM = 1 and C BM = 0. Because 77, ~ (uu + dd)/y/2 
contains two quark flavors the gluonic coefficient C BM is 
twice as large as for r] s (the normalization factor l/v2 is 
absorbed in hard kernels). 

The hard kernels T lt j(u, z) in 132|) are 



T 1J (u,z) = Y^bi(u,z)C^C%x 

i 

X (j>VM 



1 - 25 



i,3 



25 



i . 6 , 



(34) 



where Cf 1 ^ , C 1 ^ 2 are coefficients describing the content 
of final state mesons going in the n, h directions respec- 
tively, while 5ij are Kronecker deltas. The hard ker- 
nels ?2j(u, z) are obtained from the above expression 
with the replacement 1 «-> 2, while the hard kernels 
Tif (u), T2^(u) are obtained from Tu(u, z), T%j(u, z) by 
replacing bi — * Cj. For the reader's convenience, explicit 
formulas for the hard kernels are provided in Tables EJ 
El for the AS = and AS = 1 decays of B°,B- and 
B® mesons (FKS mixing scheme is used for treatment 
of 77, T}' final states, see discussion above and below Eq. 
(14111 1. Note that our phase convention for final states, 
7r + ~ ud, 7r° ~ (uu — dd)/y/2, n~ ~ du, K° ~ sd, 
K° ~ ds, K + ~ us, if - ~ sm, 7y g ~ (mm + dd)/V2, 
rj s ~ ss, differs from the one used in 0, . 

The terms in the third and fourth line of Eq. (|32|l are 

coming from Q*g) and Q^g) operators respectively. They 
contribute only if M\ is an isosinglet pseudoscalar meson, 
with hard kernels Tfj(u, z), Tf^(u) that start at NLO in 
a s (m b ) U3, (CU 



1M2 



T 3 XJ 



2 ^BM 2 



2 BM 2 C S 



(/) 



(35) 
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TABLE I: Table of hard kernels for AS = decays of B" , B° and B° into isosinglet mesons, separated respectively by horizontal 
lines. The superscript (d) in the coefficients c' d ' is not displayed. Tu and Tu are obtained with the replacement Ci — > bi. 



Mode 


± 1C 


2V 

- 1 2C 


Tr~r/ q , p~r\ q 


-7=(ci + C4) 


-Trfe — C3 + C4)+\/2(c5 — Ce) 




Cl -(- C4 




n~<j>, p~<j> 





C5 + C6 


TY~LO, p~ U! 


4=(ci + c 4 ) 


-75(02 + c 3 + c 4 ) + \/2(cs + c 6 ) 


7T^T1 f)^ 71 

'' '/siP 'Is 


__L_ { r> n IT" — /^/i ^ 

y^^C2 -+- C3 ^4J 


5 ^ 


IQ ? r /y 


1 (Co q= — C4) 


~\{C-2 — C3 + C4) — C5 + Cg 


P T 


n 

u 


— 1 4- /-v> ~\ 
^(Cs + C 6 J 


m®f , rfl, , 
7T LO,p UJ 


i(C2 T C3 — C4) 


— 2^2 + C3 + C4 ) ~ c 5 — c 6 







7)5(05 +c 6 ) 


r/qLJ, LOU! 


|(c2 T C3 + c 4 ) + c 5 =p c 6 


| (C2 + C3 + C4) + C 5 + Cfj 




2 (C2 — C3 + C4) + C 5 — C6 


3(02 — C3 + C4) + C 5 — Ce 




^(C2 - C3 + Ci) + \/2(c 5 - c 6 ) 


^(c 5 -c 6 ) 


r} s r}s 


C5 - c 6 


C5 - c 6 




{0,c 4 } 


c 5 ± c 6 


K oM c, K oM Vq 


{0, ^ C4} 


7)5(02 ± c 3 + c*) + 72(c 5 ± c 8 ) 



where Cgj^ is a coefficient describing the g/ flavor con- 
tent of M 2 , with q the flavor of spectator quark. For in- 
stance C^_ = 1, cg^ = -1/V2, while C«^_ = 0. 

The discussion of B decays into transversely polarized 
vector mesons is complicated by the presence of mb en- 
hanced electromagnetic operator, as pointed out recently 
in Ref. j54[. Namely, for the decays into neutral Vj_ 
mesons the 1/A enhanced operator 

Q^ 1 ) = rof fhtPRb^lB^^], (36) 

also contributes. Here ieB^ ^ is a purely electromag- 
netic operator related to the (n, _L) component of the 
electromagnetic field strength, defined in the same way 
as the gluonic counterpart in (|llf> . but with the replace- 
ment of QCD Wilson lines and covariant derivatives with 
the QED ones, W n — * W 7n ,iD cfl — > id^ + eA^. At lead- 
ing order in the electromagnetic coupling the operator 
(j^BTl is a result of a tree level matching of C>7 7 Q and of 

four quark operators Oi in with a photon emitted from 

a closed quark loop pS4|. 

The operator l|3(jfl leads to contributions that are mj, /A 
enhanced compared to the amplitudes for B — ► VyV|| in 
Q32JI. but which are, on the other hand, also a cm sup- 
pressed due to the exchanged photon. Numerically thus 
the contributions from (|36fl can be expected to be smaller 
than the 0(m°) terms in (|32fl . Thus at leading order the 
only contributions to B — > Vj_Vj_ can arise from nonper- 
turbative charming penguins A cc 9], possibly explain- 
ing large transversely polarized amplitudes in B — > <pK* 
decays |5^|, while the other terms are either l/rrib or 
Qn m Wf,/A suppressed (for alternative explanations see 
[22I l56j|). Therefore, a complete first order treatment 
of observables, most notably all the CP asymmetries, in 



TABLE II: Table of hard kernels for AS = 1 decays of B~ , B° 
and B° into isosinglet mesons (separated by horizontal lines 
in the table). The superscript (s) on coefficients c' s ' is not 
displayed. Tu and T2,/ follow from the replacement Ci — > bi. 



Mode 




T2( 




^(Cl+C4) 


^j(c2 T C3 + 2c 5 =F 2c 6 ) 




{ci +c 4 ,0} 


C 4 + C 5 =F C6 




75 C4 


^75(c 2 Tc 3 + 2c 5 T2c 6 ) 




{ C4 ,0} 


C 4 + C5 =F C6 







^(C2-C 3 ) 







7f( C2 + C3 ) 


VqK°,VqP° 





|(c 2 Tc 3 ) 


VqVq 


^ - ^ + cs - c 6 


a ~ 2 + c s - c e 


VsVq,<t>Vq 


-^(c 4 + c 5 =F c e) 


^5(c2-C3 + 2c 5 -2ce) 


TlsVs,<t>4> 


C 4 + C 5 =F C6 


C 4 + C 5 =F C6 


<f>Vs 


C4 + eg + C6 


C 4 + Cg — C6 


rj a Lu, <f>u> 





^(02+03+205+206) 



-B — > Vj_Vj_ decays requires an inclusion of l/rrib opera- 
tors, which is beyond the scope of the present paper and 
will not be pursued further. 

The amplitude (|32|l has the form of a convolution 
of nonperturbative light-cone wave functions 4>m(x), 
<& 9 M (x), <f) J ^(k + ) and the perturbative hard kernel and jet 
functions. With prior knowledge of the light-cone wave 
function from a fit to an unrelated experiment, a predic- 
tion for B — > MiM 2 decays can be made using a per- 
turbative expansion in as(v»%A) for the jet functions 
and in as(rnh) for the hard kernels. Alternatively, the 
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TABLE III: Table of hard kernels for AS = decays of 
B~ , B° and B® mesons respectively (separated by horizon- 
tal lines in the table) without isosinglets in the final state. 
The superscripts on c\ d ^ are not displayed for brevity. Tu 
and T\j are obtained with the replacement Ci — > bi. 



Mode 


Tic 




r 2C 


ir~iv , p~iv 


7f( Cl+C4 ) 




^c 2 — c 3 — c 4 ) 


7T p , p p 


75 ( Cl + C4 > 


1 | 


^C2 + C3 — C 4 ) 


K 0M K -M 


C4 





7T + 7V~, p + TT~ 







Cl + c 4 


K + P~, P + P~ 







Cl + c 4 


J J „0 J) 
TT IT , 7T p 


5 ( - c 2 + c 3 + c 4 ) 


l(- 


- c 2 ± c 3 + c 4 ) 


pV 


I ( — c 2 — c 3 + c 4 ) 


l(- 


- C 2 - C 3 + C 4 ) 




c 4 







■k~K+ { *\ p-K +( * ) 


Cl + C 4 









^5(c 2 TC3-c 4 ) 








nonperturbative parameters can be fit from observables 
in B — > M1M2 decays. This approach is especially useful 
at leading order in as{mb), since then the hard kernels 
Ti(,2((u) are constants, while T\j^j{u, z) are functions 
of u only. Furthermore, at this order the hard kernels 
Tfj 2J (u, z), Tft* 2 £ (u) do not contribute at all. At LO in 
as{mb) thus 

Ab^MiM 2 = 

+ fa£ BM *T K + 1^2 + AW^'^J, 
where £ BMl and 

Cf^ = JdzC^(z), (38) 

are treated as nonperturbative parameters to be fit from 
experiment. Note that in this way no perturbative expan- 
sion in as(V m bA) is needed. Eq. I|37|) has exactly the 
same form as a factorization formula obtained in Refs. 
0, ITol ] for decays into nonisosinglet mesons, but is now 
valid also for decays into isosinglets. We will use this 
form for the factorized amplitudes in the phenomenolog- 
ical analyses in Section ITTT1 

We turn now to the treatment of rj and rf states, where 
we use the FKS mixing scheme |57| . An arbitrary isospin 
zero biquark operator O can be written as a linear com- 
bination of O q ~ (uu + d(J)/v2 and O s ~ ss opera- 
tors with well defined flavor structure (here Dirac struc- 
ture is ignored for simplicity). The matrix elements of 
O = c q Oq +c s O s between 7/ ' states and the vacuum can 
then be parameterized in a completely general way by 

(0|O|7?) =c q cos (f)q(Oq) -c s sin (f> s (O s ), (39) 
(0\O\r/) = c q sin <p q (O q )+c s cos S (O S ), (40) 



TABLE IV: Table of hard kernels for AS = 1 decays of 
B~ , B° and mesons respectively (separated by horizon- 
tal lines in the table) without isosinglets in the final state. 
T2J and Tu are obtained with the replacement a — > b t . 



Mode 


Tic 




T2( 




7?( c i S) + c i S) ) 




;4 s, t4 s) ) 


gov*-, K°^p- 






7?°<*V, K°Mp° 


1 _W 






K' (t) TY+, K~ M p + 


c«+c<'> 







K -(*) K +(*) 









go(*) K a(*) 










where the four matrix elements (p\O qiS \r]^) have been 
traded for two angles <fi qtS and two reduced matrix el- 
ements (0 9jS ), all of which in principle depend on the 
structure of operator O. Phenomenologically, cf> q = cf> s = 
4> to a very good degree, with cf) = (39.3 ± 1.0)° irrespec- 
tive of the operator O [52|. In other words, if the mass 
eigenstates 77, rf are related to the flavor basis through 

r) = r) q cos <j>-T] s sin <j>, 
r) = r\ q sin <f> + r) s cos <fi, 

then to a very good approximation (i) the matrix ele- 
ments corresponding to OZI suppressed processes van- 
ish, (0|O g |7i s ) = (0|O s |?7 9 ) = 0, and (ii) the shapes of 
r)q^ s components of the wave functions do not depend on 
the mass eigenstate so that the reduced matrix elements 
(O g )„ = (O|O 9 |?7>/cos0, (O q ) v , = (O|O 9 |r/)/sin0 (and 
similarly for O s ) are independent of the final state parti- 
cle, (O q ) v = {O q ) v <. 

We will make one further assumption 

<l>ir{u,n) = <f> Vq (u,fi), (42) 

that is well respected by data. The relation 142|l is true 
for asymptotic forms of LCDA, where for p — > 00 one 
has 4>tt(u) = 0^(11) = 6uu. It can be, however, only 
approximately true for all other scales, since <f> riq {u 1 p) 
mixes with gluonic LCDA $^ (u, /j), while 0^(u,/^.) does 
not. Even so, for p as low as p = 1 GeV the relation 
(|42|l is very well respected [Hi]. The inverse moments or 
tt and r\ q LCDA for instance agree within experimental 
errors, which are at the level of 3 — 5%. 

Eq. I|42|) leads to relations between the nonperturba- 
tive function Qj) and A cc for decays into 7r,7y 9 , so that 
one can write 

C5) 9 =Cf,)+2C (J)ff , (43) 

with C(J)g new nonperturbative functions that are en- 
tirely due to contributions from interpolating gluons (the 
function £j g for instance is equal to the second term in 
(|33|l ). The decomposition (|43|l is useful only in the limit 
of exact flavor SU(3) symmetry, when 

Q*=t(J) B - (44) 
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FIG. 6: The gluonic charming penguin contributions with 
intermediate on-shell charm quarks annihilating into two 
collinear quarks going in the opposite directions, with n 
collinear quark annihilating with spectator quark and pro- 
ducing two n collinear gluons (compare also with diagrams b) 
and d) of Fig. 0). 



Similar relations for A cc will be given below. 

Using SU(3) symmetry further relations are possible. 
In the exact SU(3) limit only two £ functions are needed 
for the decays without isosinglet mesons 



C(J) 



/-Bit 



/■BK 



C 



B„K 
(•/) ■ 



(45) 



The amplitude receives contributions from SCET 
operators of higher order in l/mb, where the spectator 
quark directly attaches to the weak vertex. These higher- 
order corrections correspond to penguin annihilation in 
the diagrammatic language and do not contribute to A^) . 

At LO in 1/ rrib one further parameter is introduced for 
AS = decays into two kaons 



\KK _ aK"K 



= A 



K U K U 



(50) 



while higher order penguin annihilation contributions to 
Af? c K distinguish between the two amplitudes. 

Three additional complex parameters describe charm- 
ing penguins in AS = 1 decays of B°, B~ 



a Kit 

V2A^ r '; 



aK-tt 

1 



=.4 



Ktt 



ccg ^ 

aKj] s , AK7] a _ aK~T] s _ aK°V. 
cc,g ' cc — cc cc 



A^ = -V2A«:*° = V2A«-*° 



(51) 



where the gluonic component A^J 1 ^ has been pulled out 
for later convenience. An additional six complex param- 
eters describe charming penguin contributions in B® de- 
cays 



Furthermore, to describe all the decays into isosinglet 
mesons only the two new functions C(j)g defined in (|43[l . 
are needed. Namely, in exact SU(3) one has (cf. Eq. 



("0 



2 C(j) fl 



/-B B r} s * 



C(J)s 



(46) 



in addition to the relations , jH^l . 

Let us now discuss the nonperturbative charming pen- 
guin contributions A^[ lM2 (I37|l in the isospin limit assum- 
ing FKS mixing along with relation (|42|l . The charming 
penguins in B°,B~ decays into Trq q , Trr) s and 7T7t final 
states are parameterized in terms of four complex pa- 
rameters 



a: 



A*Vs 

cc,g 



A™ w 

A n ~Va 



-V2A£i- = y/2A%"-, (47) 



and A^2, A cc g in terms of which 



V2(A% 



cc g J 



ATT T]q — _/d 7r? ? _ 

"cc "cc "cc,gi 



AVqVq — /I" 



^"cc,gi 



(48) 



and A* c 71 =0. Here A cc g" describes the charming pen- 
guin contributions, where the n collinear quark coming 
from the annihilation of charm quarks annihilates the 
spectator quark and produces two n collinear gluons, Fig. 
El At LO in l/m b there is one additional relation 



A nr > 



(49) 



V2 



V2A2? g «(s) 
2A^(s) + 2A^(s) 
A^ K (s) + V2A^(s)^A 

A%"(8) + A%>(a) = A 



= A c 
= A 1 

— "c 



Ac 



= A c 



,K° Vq 

1 



= A 



-V2A, 

B°^K°K : 



7T K 



(52) 



where the subscript g again denotes gluonic contributions 
as before. Note that the above relations are valid to all 
orders in the as(mb) and l/mb expansions, under the 
assumptions leading to FKS mixing along with relation 

In the limit of exact SU(3) and at LO in l/mb the 
above seventeen complex parameters in (|47|) - (I52|I are re- 
lated to only two complex parameters 



A r , 



and 



Al 



Al r J = A 



Kit 



A K y 



A 



KK 



00 = A« K {s) = Al?>{s) = A£*{s), 



(53) 



Ac 



A^Vq 

cc,g 



A^Vs 
cc,g 



AKVq _ A Kr la 



A™ (s) = A™-{s) = A«f(s) = A%>-(8) 



(54) 

The same relations also apply to B decays into two vec- 
tor mesons, with the replacements r\ q — ► u>, rj s — > (f>, 
tt — > p. K —* K* , but with additional simplification since 
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the gluonic contributions i|54[l vanish. The relations ap- 
ply for each polarization of the vector mesons separately, 
with transversely polarized vector mesons receiving only 
contributions from nonperturbative charming penguins 
The relations between A cc for B — > PV decays are 
slightly more complicated, because separate nonpertur- 
bative parameters are needed if pseudoscalar or vector 
meson absorb the spectator quark. Since we will not per- 
form the phenomenological analysis of B — > PV decays 
we do not display these relations. 



C. Semileptonic B decays into isosinglet mesons 

Using the above derivations, it is fairly straightforward 
to obtain the results for vector and axial vector form 
factors in semileptonic B decays to light pseudoscalar or 
vector mesons at q 2 = 0. The V — A current (?7^(1 — lF,)b 
is matched to SCE^ LO and NLO currents 0, S IS 

HUB 



where 



Jin = — n>J '[<ln,Ui_ig$n,u 2 P Lbv\, 

_2 

3m = ^\q n ^ig^^P L b v \ 

m B 



(55) 

(56) 
(57) 

(58) 



with the ellipses denoting the remaining SCETi currents 
that either receive contributions only at NLO in as(mb) 
or contribute to the form factors at subleading order in 
1 / TTtjj. At leading order in as(ms) 



CO = Ci n = Cm = 1. 



(59) 



Since j'q, are equal to n^Q^f^d^s m fl^l 1 f° r b — * 
u, d, s transitions respectively, while jf ft is equal to 

^Qid2d2s> we can rea dily obtain the result for the 
SCETi to SCETn matching of the V-A current from the 
results obtained in the previous subsection for the four 
quark operators. The soft overlap contribution is pro- 
portional to n M , while the hard scattering contributions 
contribute equaly to terms with n M and Lorentz struc- 
ture. More precisely, the soft overlap for b — ► d decay 
equals n^Ti^d, leading to a contribution of ra B n^C, BM 
to the matrix element (Af|dy(l - js)b\B) at q 2 = 0, 
with similar results for the hard scattering contributions. 
Using a definition of form factors 



(p\rf(i 



l5 )b\B) =C BP [m B n»f BP (0) 



m B 



^(/f p (o) + / BP (o))], 



(60) 



where the relations p B = m B v^ , and p P = m B n^ jl have 
been used, this gives for the form factors of B — > P tran- 
sition at maximal recoil at LO in l/m& and as(mb) 



f? P (0) + f BP (0) 



BP i 



(61) 
(62) 



Quite similarly one obtains for the form factors at q 2 = 
in B — > V transition 



m B 



MO) = 7T^(M0) - 4,(0)) = A 3 (0) = C BV + Cj 



2my 



and 



A 2 (0) + 2m v m B 



d{A 3 - A 



dq 2 



m B 



(63) 



(64) 



where the standard definition of form factors has been 
used (and can be found e.g. in .58]). The coefficients 
C BP in the definition of B ~ » P form factors (jfif)|> and 
equivalent coefficients C BV in the definition of B — ► V 
form factors, take care of the flavor content of the fi- 
nal state M. For instance for rj q ,uj these are C B M — 
C B ° M = l/y/2, while for they are C B ~ M = 

C = 1. The expressions for (j parameters in terms of 
jet functions are given in Eq. O. Note that the derived 
relations between form factors and the SCET nonpertur- 
bative functions are valid also for isosinglet mesons and 
in this sense extend the previous discussions of B — > P, V 
form factors [3^, E3, EH |43, |4^, EH E^, E3| ■ In partic- 
ular, the B — > r)q form factors /±(0) receive the gluonic 
contribution as can be seen from 

III. PHENOMENOLOGY 

We now apply the LO (in 1/to& and as(rrib) expan- 
sions) factorization formula Eq. I|37l) to B decays into 
pseudoscalar mesons. 



For the nonperturbative parameters A 



Ah M-2 



and ( 



BMi 
(J) 



exact SU(3) relations JSJ-gBJ, J33J|, JSU will be used, 
leading to four independent real parameters ((j), C(J)g 
and two complex parameters A cc , A ccg describing non- 
perturbative charming penguins. The parameters C(J)g 
and A ccg correspond to gluonic contributions and are spe- 
cific to the decays into isosinglet final states. We will thus 
first determine Gj) and A cc using a x 2 " nt to observables 
in B — ► 7T7T, B — ► ttK decays and then fix the remain- 
ing parameters C(J)g an( i A ccg from a separate x 2_n t to 
observables in B — ► rf- V, B — ► rf^K decays. 

The data used in the analysis are from HFAG, Summer 
2005 compilation [5{j, apart from the observables where 
different experiments do not agree, in which case the er- 
rors are inflated according to the PDG prescription [6(| . 
The observables used throughout the analysis are the CP 
averaged decay widths 



W\ i, 

8tttti b 2 



(65) 
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with \p | = 771^/2 + 0(m 2 M /m 2 B ) the three momentum of 
light mesons in the B rest frame, and an additional factor 
of 1/2 on the right hand side, if two identical mesons 
are in the final state, while the abbreviations are Af — 
Ag^j, Af — A B ^j. Then the direct CP asymmetries 



Af 



\Af\ 2 -\A f \< 
\Af? + \Atf 



(66) 



and the additional observables, Sf and Hf, that can 
be extracted from time dependent decays of neutral B 
mesons are (restricting / to have definite CP) 



T(B° q (t) - /) = e- rt T(B q - /) 



cosh 



AYt\ 



Hf sinh \— —J - _4y F cos(Amf) - 5/ sin(Ami) 



(67) 



where Am = tuh — tul > 0, T is the average decay 
width and Ar = Th — ^l the difference of decay widths 
for heavier and lighter B® mass eigenstates. The time de- 
pendent decay width T(B®(t) — > /) is obtained from the 
above expression by flipping the signs of the cos(Amt) 
and sin(Amt) terms. 

Since (Ar/r)s d <C 1 in the B® system, only measure- 
ments of the parameter 



Im 



St 



e-^Af(A f Y 



l^/l 2 + MM 



(68) 



are experimentally feasible in the foreseeable future. 
In the above relation (|68J) the phase convention with 
argV c fc = argils = arg(-V cd ) = 0, so that 8 = 
axg(-V* b V t * d V cd Vtb) = a,rg(V* d V tb ) was employed 

In the B s system, on the other hand, we expect a much 
larger decay width difference (AT/T)b s = -0.12 ± 0.05 
within the Standard Model [62], while experimentally 
(AF/F) Bs = -0.33t°;?? H3, so that both 



(S f ) Bs = 2- 



Im e +i2e A f (A f y 



l^/l 2 + l^/l 2 



(69) 



and 



{Hf) Bs = 2- 



Re 



Af(A f Y 



(70) 



might be experimentally accessible |63| . Thus predictions 
for both innj and (|7UJ) will be given in Section [TTTdI In 
writing (|69|l and 1|70|) the same phase convention as in 
(P| was used, with e = axg(-V cb V ts V* 8 V* b ) jEJEl- 

In our numerical estimates we use preferred values 
of Ref. |H3 for inverse moments of LCDA at 1 GeV 
(x- 1 )^ = 3.3 and {x~ x ) K = 2.79, (x~ x ) K = ((1 - 
Xs^^k = 3.81, and take (z" 1 )^ = (z -1 )^, = (x^ 1 )*. 



Note that this choice respects the sum rule for inverse 
moments 661 



{x~\ + Six- 1 ),, = 2((x- 1 ) K + {x- x ) K ). 



(71) 



The inverse moment of the pion LCDA also agrees with 
the experimental determination (x' 1 )^ = 2.91 ±0.54 ob- 
tained from CLEO data on the pion's electromagnetic 
form factor in Ref. [§3 . The values of the CKM el- 
ements and the CKM unitarity triangle angles 7 and 8 
are taken from CKM fitter, Summer 2005 update • In 
particular, \V ub \ = (3.899±0.10) x 10~ 3 , 7 = 58.6°±6.4° 
and 8 = 23.22° ± 0.75°, while for the weak phase in the 
B° - B° mixing e = 1.04° ± 0.07° 68]. These values 
agree with the ones obtained by the UTFit collaboration 
|69| . Note that the constraints on CKM angle 7 (a) are 
obtained using an isospin decom pos ition of time depen- 
dent B — > 7T7T, pir, pp decay s l70l |7l| (for a discussion of 
isospin violating effects see |72j ) and from B — > DK de- 
cays E3 • No use of theoretical inputs from a l/m b 
expansion has been made at this point (such a possi- 
bility for using SCET to facilitate extraction of a from 
B — > 7T7T has been discussed in jz3). Also, B — > mr data 
are not restrictive at present in the determination of a, 
so essentially no use of B — > PP data has been made to 
fix the above CKM parameters to be used in our anal- 
ysis. For the decay constants we use f„ = 131 MeV, 
f B = 218 ± 23 MeV jz3 and f Vq = 140 ± 3 MeV, 
f Vs = 176 ±8 MeV 

Since our phase convention for the ir states, 7r + ~ ud, 
7T° ~ l/y/2(uu — dd), ir~ ~ du and the kaon states K° ~ 
ds, K° ~ sd, K + ~ su, K~ ~ us differs from the one 
used in Refs. [tj 0] the corresponding hard kernels are 
gathered in Tables IITlI and ITVl 



A. Analysis of B — *• nn and B — > nK decays 

Let us start the analysis with a discussion of the 
AS = B — > 7T7T and AS 1 = l/3-> 7T.fr decays. Our anal- 
ysis differs in the treatment of errors and in the way the 
SCET nonperturbative parameters C(J)i A cc are deter- 
mined from LO SCET analysis of the same decay modes 
111 Refs. Here a combined X 2 ~fit to available ex- 

perimental information will be made to determine the 
SCET parameters, while in ^(j only a subset of modes 
was used for this purpose. The values of the SCET pa- 
rameters that will be determined in this section will then 
be used in subsequent sections; in the discussion of de- 
cays into final states with 77' ' mesons in Section IIII Bl 
for predictions on S parameters in penguin dominated 
modes in Section lill CI and for prediction of observables 
in B® decays in Section fill Dl 

The majority of the differences between B — > mr and 
B — * ttK decays can be explained by a CKM hierarchy of 
different contributions. Here and in the rest of the analy- 
sis we will split the amplitudes into "tree" and "penguin" 
contributions according to the CKM elements. Using the 
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unitarity relation A 



(d,s) 



-A 



A« we define 



Ab^ } = ^u' s) T§^ f + I'n ... (72) 

for any process B — > /. In addition we will also use 
the nomenclature where tree contributions (without quo- 
tation marks) will denote insertions of 0\ 2 @ opera- 
tors and will be present only in "tree" amplitudes Ts—ff, 
charming penguin contributions will denote insertions of 
0\ 2 and will be present only in "penguin" amplitudes 
Ps^f, while QCD penguin and EWP contributions will 
denote insertions of the operators 03,...,6 and O7 10 re- 
spectively and will contribute to both Ts^f and Ps-tf- 
In AS = 1 decays, such as B — > wK, there is a strong 
CKM hierarchy between the two terms in (|72[1 since 

|Al s) | ~ 0.02| A£ s) I . No such CKM hierarchy is present 
in AS — decays, where both terms are of the same or- 
der in the Wolfenstein expansion |Al d) | ~ \\i d) \ ~ A 3 [77|. 
One therefore expects sizeable direct CP asymmetries in 
AS* = processes and much smaller ones in AS = 1 
decays. 

There are a number of other observations that can be 
made using the LO SCET expression (|3T|l even before 
the nonpcrturbative parameters of SCET are determined 
from the data. Because of the intriguing discrepancies 
between theoretical expectations and experiment in the 
AS* = 1 transition, we will focus in this section mainly 
on B — > ttK decays. We start with the isospin decompo- 
sition of B — > nK decay amplitudes [73, |ZS H3 



Ab°^tv+K ^1/2 — ^3/2 — A{, 2 , 



1/2 ^3/2 ^1/2 
1 .n i= 1 



(73) 



Ab-^>k- =^A° 1/2 - V2Al /2 + ^=A\ /2 , (74) 
Aso^tvORo = - "Tg^i/a - ^^3/2 + "y^ 1 / 2 ' ( 75 ) 



Ab-^tx- K° —A1/2 + A\/ 2 + A i / 2 , 



(76) 



with / and AI in the reduced matrix elements Af 1 de- 
noting the isospin / = 3/2,1/2 of the final states and 
AI = 1,0 denoting the isospin content of the weak 
Hamiltonian. If electroweak penguin contributions are 
neglected, the reduced matrix elements A^ 2 and A\ , 2 
receive contributions only from tree operators 0\ 2l while 
A\/2 rece i ves both penguin and tree contributions. Ex- 
plicitly, at LO in SCET 



A 1 

^3/2 



3 V2^ 1 

uk BK (c^ 



K 



*Btt 



.(«) 
-1 



)]}■ 



(79) 



(s) (s) 

where SCET Wilson coefficients c\ ,b\ are to be under- 
stood as already convoluted with LCDA. This amounts 
to a replacement 0^2,3 — > ±m,B / (xj^) k , in b[ s ^ that are 
multiplied by fx and a replacement u>2,3 —> i7ns/(a; 

(s) 

in 6 2 that are multiplied by f n . 

The dominant term in B — > J^tt decays is the charm- 
ing penguin term A^ 7r , which is A^/A« enhanced over 
tree contributions. Since it arises from insertions of 0\ 2 
operators it is also larger than QCD penguins by a fac- 
tor Cias(2m c )/ max(C3, . . . , Cq) ~ 10. The charming 
penguin contribution has AI — and is thus present 
only in A®, 2 l)77p. Similarly, the reduced matrix elements 

A3 , 2 1 ^ do not receive contributions from QCD penguins 
since these are AI = operators. Note also, that the 
presence of the inverse moment (x ) n g ~ 3 in hi lifts 
the color suppression of tree operators in B° — > K°ir°. 
Similarly the color suppression of EWP contributions to 
B° — > 7r + X _ and B~ — > 7r _ X° is lifted as well. Isospin 
decomposition also leads to the relation 



A, 



1 



(80) 



+ ^S°-i-if 7r ! 

valid to all orders in 1/rrih and as [x3jE3.IH3| • As pointed 
out recently in Ref. [8IJ this relation furthermore receives 
corrections from isospin breaking which are of only sec- 
ond order numerically, with corrections to penguins that 
are linear in m u d/A canceling exactly. 

The dominance of the A^? term in B — ► Kir ampli- 
tudes leads to the approximate relation between branch- 
ing ratios 



Bi 



-rr+K- 



Br 



TT-K° 



2Br 



TT°K- 



2Bvx», (81) 



that is well obeyed by the data. The corrections to 
these relations come from A^ 2 3 ^ 2 reduced matrix el- 
ements (|78f) . <|79[) that receive only contributions from 
\Xu^\ ~ 0.02| a£ s ^ I suppressed tree operators or from 
EWP. To study them it is useful to construct ratios of CP 
averaged decay widths in which the dependence on A^ 
cancels to first approximation. Following the notation in 
the literature we define 85 



A? /2 =mi^§{AW^ + - 



f_K 

2 



+ C f (6«+2&< s >)]}, 
! m 2 B G F ( 

+ 2U\C BK (c¥-ci s) )+Cf K (b 



(77) 



BK( h (s) _ 



>]}■ 



(78) 



R = W^£) B - aM±am - (82) 
»<=> Z- im±om - (83) 

_ l T(B°->K-n+) Exp 
Rn - 2 f (50 ^ W) " °- 82±0 - 08 ' ( 84 ) 

where the experimental information on the ratios has also 
been displayed. For the reader's convenience we will also 
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TABLE V: Predicted CP averaged branching ratios (xlO -6 , 
first row) and direct CP asymmetries (second row in each 
mode) for AS = and AS = 1 B decays (separated by hori- 
zontal line) to two nonisosinglet pseudoscalar mesons. The er- 
rors on the predictions are estimates of SU(3) breaking, 1/mj 
corrections and due to errors on SCET parameters, respec- 
tively. 



Mode 


Exp 


Theory 


B° 




5.0 ±0.4 


5.4 ± 1.3 ± 1.4 ±0.4 






0.37 ± 0.23 a 


0.20 ±0.17 ±0.19 ±0.05 


B° 


- ttV 


1.45 ± 0.52 b 


0.84 ± 0.29 ± 0.30 ± 0.19 






0.28 ± 0.40 


-0.58 ±0.39 ±0.39 ±0.13 


B~ 


- 7r°7r- 


5.5 ±0.6 


5.2 ± 1.6 ± 2.1 ± 0.6 






0.01 ± 0.06 


< 0.04 


B~ 


-> K°K~ 


1.2 ± 0.3 


1.1 ±0.4± 1.4 ±0.03 


B° 




0.96 ±0.25 


1.0 ±0.4 ±1.4 ±0.03 


B° 


— > 7T K 


11.5 ± 1.0 


9.4 ±3.6 ±0.2 ±0.3 






0.02 ±0.13 


0.05 ±0.04 ±0.04 ±0.01 


B° 


-> K--K+ 


18.9 ±0.7 


20.1 ±7.4 ±1.3 ±0.6 






-0.115 ±0.018 


-0.06 ±0.05 ±0.06 ±0.02 


B~ 


->K-TT° 


12.1 ± 0.8 


11.3 ± 4.1 ± 1.0 ±0.3 






0.04 ± 0.04 


-0.11 ±0.09 ±0.11 ±0.02 


B~ 




24.1 ± 1.8 c 


20.8 ±7.9 ±0.6 ±0.7 






-0.02 ± 0.05 d 


< 0.05 



PDG scaled errors a (S = 2.3), b (S = 1.8), C {S = 1.4), 
d (S = 1.5). 



discuss the ratio Rqo = R/R n proposed in 01 

The deviations from 1 are experimentally only at the 
level of at most 2a. Also, of the CP asymmetries only 
A^L^ + is well measured, so that the values of the SCET 
parameters (PJf , and the phase of A^r cannot at 
present be reliably extracted from B — > Kix experi- 
mental data alone. We thus impose SU(3) symmetry 
A cc = AH = Kc, C(J) = Cfjf = Cfj) and construct X 2 
from observables in B — > 7T7t and B — > irK decays. 

In B — > 7T7T decays there is experimental information 
on seven observables: the time dependent CP asymme- 
try S v +„-, the three CP averaged decay widths T(B — > 
7r+7r-),f(B -> 7r°7r°),f(B -> 7r~7r°) and three direct 
CP asymmetries A°X V - , ■A^o ^ P ,o- This latter is not 
used in the fit since B~ — > tt~tt° is a AI = 3/2 pro- 
cess and thus does not receive QCD or charming penguin 
contributions, so that strong phases are generated only 
at NLO in as(m{,), while at LO the asymmetry is zero 
irrespective of the SCET parameters. 

In addition, the following observables in B — > Ktt 



TABLE VI: Predictions for the CP violating S parameters. 
The errors on the predictions are estimates of SU(3) break- 
ing, 1/rnt, corrections and errors due to SCET parameters, 
respectively. 





Mode 


Exp 


Theory 


B° 


— » 7T - 7T + 


-0.50 ±0.19° 


-0.86 ±0.07 ±0.07 ±0.02 


B° 


- ttV 




0.71 ±0.34 ±0.33 ±0.10 


B° 




0.31 ±0.26 


0.80 ±0.02 ±0.02 ±0.01 



Error scaled according to PDG (S=1.5). 



decays are used in the x 2_n t ; the four CP averaged 
decay widths T(B K°ir°), f (B -> K-tt + ),T(B -> 
K-tt°),T(B -> K°ir~), but only three direct CP asym- 
metries A^^otA^^+tA^^o- The prediction for the 
remaining direct CP asymmetry AfE, _ can receive large 
corrections at NLO in as(nib) from terms of the form 
Al s 'Ci j 2Q:s( TO 6)- These can be comparable in size to LO 

(s) 

terms proportional to A« which come entirely from QCD 
penguin operators. Also, the experimental information 
on S Ks7r o is not used in the x 2_ nt, and will be discussed 
separately in Section IlII CI 

From the x 2_n t to the B — + itit, Ktt data we then ob- 
tain 

C = (7.3 ± l. 8 )*10-* 
C.7 = (10.3 ± 1.6) x 1CT 2 , 

and 

\A CC \ = (46.8 ±0.8) x 1(T 4 GeV, 
arg(A cc ) = 156° ±6°, ''^ 

with x 2 /d.o.f. = 44.6/(13 — 4), where the largest discrep- 
ancies are in A^o K - , A^- J(+ as can be seen from TablelVl 
This very high value of \ predominantly reflects the fact 
that the expected theory errors coming from NLO l/m& 
and as{mb) terms and from SU(3) breaking are larger 
than experimental errors. If the estimates for these er- 
rors, to be discussed below, that are given as first and 
second errors on the theoretical values in Table ]V\ are 
added quadratically to experimental errors in the defini- 
tion of x 2 , the resulting value is x 2 /d.o.f. = 8.9/(13 — 4) 
(rVd.o.f. = 15.3/(13 - 4) if SU(3) breaking errors are 
taken to be correlated). The extracted values of SCET 
parameters (18611 , (|87|l agree within errors with similar ex- 
tractions from only 7T7t data or a combination of 7T7T and 
ttK data without modes that depend on that were 
performed in Ref. [Toj . 

Using the above values for the SCET parameters one 
can predict CP averaged decay widths and direct CP 
asymmetries with the results listed in Tables |V] and 
IVII These results were obtained in the limit of ex- 
act SU(3) and as such an error of 20% is introduced 
as an estimate of SU(3) breaking effects in the relation 
C(j) = (P/S = (Pn- Similarly an additional 20% error on 
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the magnitude and 20° error on the strong phase is intro- 
duced due to SU(3) breaking and 1/m^ corrections in the 
relation A cc = = . These variations result in 

the first theoretical error estimate in Tables IVII The 
second error is an estimate of the remaining 1/mb and 
ots{mb) corrections which we take to have a magnitude 
of 20% of the leading order contributions proportional 
( f) ( f) 

to A« and AJ CKM elements and assume that they 
introduce an error of 20° on the strong phase (with the 
exception of A cc in AS = decays, where the error of 
5% on the magnitude is assigned due to insertions of 
03,...,io operators leading to nonperturbative charm con- 
tributions proportional to xlf^ instead of X c d \ while in 
AS = 1 decays this effect is negligible). 

A different approach is used for the observables in 
which the subleading corrections are expected to be 
anomalously large. This can happen in the decay modes 
in which the leading contributions are not proportional to 



„(/) u(f) 



, b\ J ' or c 2 , b 2 , so that l/m& or as(m,b) suppressed 
contributions may lead to 0(1) corrections. For instance, 



a Gf 2 

A B-^K^- = -J=™ 



b{/k 



-X ( ^A 



ttK\ 
cc J ' 



where the meaning of SCETi Wilson coefficients c\ , b\ 

is the same as in Eqs. The A « } P art 01 

the B~ — > K°tt~ can receive corrections of the form 
Ci.2(Xs( m b) which can compete in size with A«^ part of 

( s) { s) 

c\ , b\ in LO hard kernel. Thus only a conservative 
upper bound on the CP asymmetry is given in Table [V] 
Similar reasoning holds for B~ — » K a K~ decay 



.4, 



>K°K~ 



Gf 

: 7f 



f-BK Id) r BK h {d) 



A^Af/}, 



(89) 



and for B° — > K°K° decays, whose amplitude is equal 
to (|S9|) at LO in 1/mfe and as{rrib). Again, for the 
"tree" part of the amplitude we can expect large NLO 
corrections from Ci t 2Cts( m b) terms. Because there is 
no CKM hierarchy between "tree" and "penguin" am- 
plitudes the predictions for these two modes are even 
more uncertain, so that we do not give any bound on 
CP asymmetry, while the estimate of as(rrib) correc- 
tions to branching ratios is taken to be the same as for 

5° -> TT-TT+. 

Many of the errors cancel to a large extent in the ra- 
tios of the decay widths, which are then predicted more 
precisely than the individual rates 



R- 


lT h. 


(3.7 ±1.5 ±3.9 ±2.1) x 1(T 2 , 


(90) 


R c - 


ll' 


(8.8 ±2.3 ±6.9 ±1.2) x 1CT 2 , 


(91) 


Rn — 


ll' 


(6.9 ±2.0 ±7.5 ±0.8) x 10~ 2 , 


(92) 


Roo — 


lT h. 


(-3.0 ± 0.9 ± 3.2 ± 1.4) x 10~ 2 , 


(93) 



with the errors estimating the SU(3) breaking, l/m,b and 
as (nib) corrections, and the errors due to uncertainties 
on SCET parameters respectively. Note that even though 
the SCET parameters were determined using these data 
as well, the agreement between predicted ratios and the 
experimental values (|82(l - (|85|l is far from impressive. Ex- 
perimentally i?, R n < 1 with about a 2a difference from 
the above expectations. 

An important input to the above theoretical predic- 
tions was provided by mr data using SU(3) symmetry. If 
instead no itit data is used and the SCET parameters are 
determined solely from the B — > irK decays, the theo- 
retical and experimental values of the R ratios l|82|l - l|84() 
would agree within experimental errors, but with values 
of SCET parameters that differ significantly from Q86p. 
with C a factor of 4 larger, while (j even flips sign! 2 This 
leads us to two conclusions, that (i) the SCET expan- 
sion on itself is not in contradiction with ttK data and 
(ii) without permitting extremely large SU(3) breakings 
that allow even for a change of sign for the nonperturba- 
tive parameters, there is a discrepancy between data and 
theoretical expectations. An independent check on the 
validity of the SCET 1/mf, expansion can be provided 
with a phenomenological analysis of B — > ttK data using 
diagrammatic decomposition, where no SU(3) is assumed 
but annihilation topologies are neglected [68ll83| . 

Especially interesting is the difference between R n and 
Rc pi. l85| . Defining the "tree" and "penguin" con- 
tents of Aj analogously to the general decomposition 
Eq. G21 



Af 1 



Xi s ~>Tj AI + AWP/ 



A/ 



(94) 



and using the fact that dominates the amplitudes 

SO that Ac ^ u "^3/2 1/2' 

A^T^ and A^P^ » 
A^Pg 1 ^ jy 2 (this latter hierarchy follows from the fact 
that Py 2 receive only EWP contributions and are 
thus smaller than charming penguin contributions in P^ 2 
that arise from insertions of 0\ 2 ), then completely model 
independently to first order in small parameters 



Rc — Rn 



The higher order corrections to the above relation 
amount to a difference 




Th 



(R c - R n ) =' (1.8 ± 0.9 ± 0.9 ± 0.4) x 10" 



(96) 



in the theoretical predictions l|83|l . I|84|) . This is to be 
contrasted with the difference between experimental val- 
ues (R c — Ru)ex P . = 0.18±0.13. If this difference persists 



2 The result of this fit is f = 
Acc = (48.8 ± 1.1) x 1CT 2 , 
unacceptably small Br(B~ - 



0.26 ± 0.05, Cj = -0.19 ± 0.05, 
arg Acc = 131 ± 13°, which gives 
^n~n°) = (0.4+° 9 4 ) X 10- 6 . 



16 



with reduced experimental errors, it will be very difficult 
to explain in the Standard Model. A possible explana- 
tion due to isospin violating new physics that changes 
EWP contributions and thus enhances the Ph 2 term 
in l|95l) has been extensively discussed in the literature 

That the ttK data alone are not inconsistent with the 
theory is well demonstrated for instance by the sum of 
CP averaged decay widths 



2f (B° -» K°n°) - T(B° -> K-n+)+ 



(97) 



first discussed by Lipkin [8|j and by Gronau and Rosncr 
|9C| . The sum (|9T)l does not depend on A\ , 2 and thus on 

A*? . In terms of the R ratios the sum 1|97(1 is 



the corrections are subleading in the l/m^ expansion 
|93| . Using isospin decomposition (|77Jl - i|7t?|) and defining 
"tree" and "penguin" terms of the corresponding reduced 
matrix elements <|94ll one has 

A E = -24 Im [A«A«*] Im [{Tfo-T^Pfo+T^Pfa] 

(104) 

The "penguin" terms Ph 2 ]/ 2 receive only EWP contri- 
butions, while the "tree" terms Th 2 l , 2 are a sum of tree 
and EWP contributions, with tree contributions domi- 
nating due to larger Wilson coefficients. At leading or- 
der in l/m,b and as(rrib) expansion the strong phase is 
nonzero only due to A^ c w . Thus at this order Th 2 

and P\j 2 1/2 are rea l ( m our phase convention) so that 
A^ = to the order that we are working. 



AL = R 



oo 



R + R c -1. 



Using the expansion to first order in small parameters 
/ P i% T 3/2 i/2/ P i°/2 as in © for i?oo and R 



P 1 

3/2,1/2 



R 



00 



R=l-4 




(99) 



(100) 



it is easy to check that AL is only of second order in small 
parameters in accordance with the fact that all interfer- 
ence terms with A®, 2 in (|97|l cancel. In AL the depen- 
dence on A^ drops out completely leading to the value 

of AL that is \\ { u )2 1 \{ s]2 )\ ~ A 4 = 2 • 10" 3 CKM sup- 
pressed. Experiment at present is consistent with van- 
ishing AL at one a 



Ai E = p ' 0.19 ±0.14, 



(101) 



while the theoretical prediction using LO SCET expres- 
sions and JHSJl, l|57|l. is 



Th 



AL = (2.0 ± 0.9 ± 0.7 ± 0.4) x 10" 



(102) 



Another very precisely predictable quantity is the sum 
of partial decay differences Ar = T(B -> /) - T(B -> f) 



A E =2AT(B- -► K-ir°) - AT(B~ - 
+ 2Ar(B° -> K"tt°) - AT{B° 



K°TT-) 



(103) 



In the hmit of exact isospin and no EWP A^ vanishes 
[9ll I92I |9^|. However, even in the presence of EWP, 



B. Decays into isosinglet states 

At present the experimental data on the decays with 
rfi > are not abundant. Of the 55 observables describing 
the complete set of B° , B~ and B® decays to two body 
final states with rp \ only 11 have been measured so far. 
We will thus assume the SU(3) relations lO-lO, 
(|54|l between SCET parameters as in the previous subsec- 
tion and furthermore use the determination of the SCET 
parameters Q nand A cc from the tttt and ttK data given 
in Eqs. I|86l) . 1)8 7|) . The remaining parameters specific to 
isosinglet modes, C(J)gi describing gluonic contributions 

to B — » r}(' form factors, and A ccg , describing the glu- 
onic parts of charming penguin, Fig. are then fixed 
from a y^-fA to observables in the isosinglet modes. 

As in the previous subsection only CP averaged decay 
widths and direct CP asymmetries are used in this de- 
termination, while the discussion of Sk s v' * s relegated to 
Section Till CI This leaves four observables from AS* = 
decays: two CP averaged decay widths and two direct 
CP asymmetries in B~ — » 7r - ?/ and six observables 
in AS = 1 decays: three CP averaged decay widths 
and three direct CP asymmetries in B° — > K°r]' and 
B~ — > K~r)( ' modes. 

In the above modes the functions £ g and £j g enter in 

the combinations Cs+&i Os an d c^Q+b^\j g , with 
the latter being numerically much smaller (cf. JSJl, i|15fl . 



(|16|I % 1. Since c\ ~ b\ J ' we expect the combination Q T 
Cj 9 to be relatively well determined from the data, with 
the orthogonal combination only poorly constrained. We 
thus define 



,(/) 



C± = Cg ± Off, (105) 

which are then fit from the data. We find two sets of 
SCET parameters that minimize y 2 : 
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TABLE VII: Predicted CP averaged branching ratios (xlO -6 , first row) and direct CP asymmetries (second row for each mode) 
for AS — and AS — IB decays (separated by horizontal line) to isosinglet pseudoscalar mesons. The Theory I and Theory II 
columns give predictions corresponding to Solution I, II sets of SCET parameters. The errors on the predictions are estimates 
of SU(3) breaking, l/m^ corrections and errors due to SCET parameters, respectively. No prediction on CP asymmetries is 
given, if [—1, 1] range is allowed at la. 



Mode 


Exp. 




Theory I 


Theory II 


id 


— > TV Tj 


4.3 ±0.5 (S = 


1.3) 


4.9 ± 1.7 ± 1.0 ±0.5 


5.0 ± 1.7 ± 1.2 ±0.4 






-0.11 ±0.08 




0.05 ±0.19 ±0.21 ±0.05 


0.37 ±0.19 ±0.21 ±0.05 


ID 


1 

— > TV Tj 


2.53 ±0.79 (5 


= 1.5) 


2.4 ± 1.2 ±0.2 ±0.4 


2.8 ± 1.2 ±0.3 ±0.3 






0.14 ±0.15 




0.21 ±0.12 ±0.10 ±0.14 


0.02 ±0.10 ±0.04 ±0.15 


B° 


-> Tr°n 
> ty rj 


< 2.5 




0.88 ±0.54 ±0.06 ±0.42 


0.68 ±0.46 ±0.03 ±0.41 






— 




0.03 ±0.10 ±0.12 ±0.05 


-0.07 ±0.16 ±0.04 ±0.90 


B° 


-> Tr°n' 


< 3.7 




2.3 ±0.8 ±0.3 ±2.7 


1.3 ±0.5 ±0.1 ±0.3 






— 




-0.24 ±0.10 ±0.19 ±0.24 


— 


ID 


— > 7777 


< 2.0 




0.69 ±0.38 ±0.13 ±0.58 


1.0 ±0.4 ±0.3 ±1.4 






— 




-0.09 ±0.24 ±0.21 ±0.04 


0.48 ±0.22 ±0.20 ±0.13 


ID 


— > r\r\ 


< 4.6 




1.0 ± 0.5 ± 0.1 ± 1.5 


2.2 ±0.7 ±0.6 ±5.4 






— 




— 


0.70 ±0.13 ±0.20 ±0.04 


B° 


-> n'v' 


< 10 




0.57 ±0.23 ±0.03 ±0.69 


1.2 ±0.4 ±0.3 ±3.7 












0.60 ±0.11 ±0.22 ±0.29 


B° 


-> K°r)' 


63.2 ±4.9 (S = 


= 1.5) 


63.2 ±24.7 ±4.2 ±8.1 


62.2 ±23.7 ±5.5 ±7.2 






0.07 ±0.10 (S 


= 1.5) 


0.011 ± 0.006 ± 0.012 ± 0.002 


-0.027 ± 0.007 ± 0.008 ± 0.005 


B° 


— > K°rj 


< 1.9 




2.4 ±4.4 ±0.2 ±0.3 


2. 3 ±4. 4 ±0.2 ±0.5 










0.21 ±0.20 ±0.04 ±0.03 


-0.18 ±0.22 ±0.06 ±0.04 


B~ 


-> K- v ' 


69.4 ±2.7 




69.5 ±27.0 ±4.3 ±7.7 


69.3 ±26.0 ±7.1 ±6.3 






0.031 ±0.021 




-0.010 ± 0.006 ± 0.007 ± 0.005 


0.007 ± 0.005 ± 0.002 ± 0.009 


B~ 


-> K- v 


2.5 ±0.3 




2.7 ±4.8 ±0.4 ±0.3 


2.3 ±4.5 ±0.4 ±0.3 






-0.33 ±0.17 (S = 1.4) 


0.33 ±0.30 ±0.07 ±0.03 


-0.33 ±0.39 ±0.10 ±0.04 



Solution I: 



arg(A Cl 



(-9.9 ±2.4) x 1CT 2 , (106) 

(-3.5 ± 14.6) x 10" 2 , (107) 

(35.8 ± 1.9) x 10" 4 GeV, (108) 

-109° ±3°, (109) 



with x 2 /d.o.f. = 25.0/(10 — 4), where the largest discrep- 
ancy with data is in «4^L . The value of \ 2 is reduced 
to x 2 /d.o.f. = 7.6/(10 — 4), if theoretical errors due to 
SU(3) breaking and estimated NLO corrections are added 
quadratically to experimental errors in the definition of 
X 2 - The second solution, on the other hand, 
Solution II: 

C+ = (-6.6 ±4.3) x 10~ 2 , (110) 
(- = (-11.2 ± 28.7) x 10" 2 , (111) 
\A ccg \ = (36.2 ± 2.2) x 10" 4 GeV, (112) 
arg(v4 CC5 ) = 68 ±4°, (113) 

hasx 2 /d.o.f. = 40.8/(10-4) or x 2 /d.o.f. = 5.4/(10-4), 
if theoretical errors are added in the definition of % 2 . The 
largest discrepancies with experimental data in this case 



is in Az^L while the prediction for AE^- agrees well with 
data in contrast to Solution I. 

The strong phases of the gluonic charming penguin 
in the two solutions lie in opposite quadrants, while the 



values of \A 



CCt g 1 and agree between the two solutions. 



The gluonic contribution to the B — > rj^> form factors, 
Cg + Os) is similar in size to £ and £j in l|86l) as expected 
from SCET counting, Using Eq. (glj we find in the SU(3) 
limit and at LO in l/m& and as (to?,) 



f+ nq (0) 



(-2.3 ±4.8) x 10~ 2 , 
^ (4.5 ±8.6) x 10~ 2 , 

'(-9-9 ±2.4) x 10 
(-6.6 ±4.3) x 10 



-2 
-2 



(114) 
(115) 



to be compared with /f "(0) = 0.176 ± 0.007, that is ob- 
tained using the results of tttt, ttK fit JSSJ. The upper 
(lower) rows in l|114|) , (|115fl correspond to values in Solu- 
tion I (Solution II), where only experimental errors due 
to the extracted SCET parameters are shown. Because of 
the large experimental uncertainties, the gluonic contri- 
butions to the form factors are still consistent with zero 
at a little above the la level in Solution II. The gluonic 
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charming penguin A ccg on the other hand is shown to be largest contributions to the amplitudes are thus 
nonzero in both sets of solutions and is of similar size to 
A cc in (|87() in agreement with SCET counting. 



That the gluonic contribution is of leading order in 
X/mb has already been recognized in the context of QCD 
factorization. In the phenomenological analysis of Ref. 
[lflj . the gluonic contributions to the B — > rj form fac- 
tor are proportional to Fi , a parameter not known from 
other sources and given the rather arbitrary values of 
F-2 = 0,0.1. With these values, the gluonic contribution 
accounts for 0%,40% of the B — > rj form factor with 
constructive interference between the gluonic and the re- 
maining contributions. This can be compared with our 
analysis where, in the B — > rj form factor, destructive in- 
terference between C(J)g and C(j) terms is found with the 
gluonic contribution from Q g + (j g 2.1 (1.4) times larger 
than the contribution from C + O m Solution 1(H)). 

The predicted branching ratios and CP asymmetries 
using the above values are compiled in Tables IVIII and 
IVIIII The errors due to SU(3) breaking and l/rrif, or 
as(r/j(,) corrections are estimated in the same way as in 
previous subsection. An error of 20% and a variation on 
charming penguin strong phase of 20° is assigned to re- 
lations 1)43(1 - 1(46(1 and ((531) . ((54(1 giving the first error esti- 
mate in the Table ED The remaining l/m^ and as(mb) 
errors, listed as second error estimates in Table Ivnl are 
obtained by varying the size and strong phase of leading 



order amplitudes proportional to A« or XY' by 20% and 
20° respectively. 

A prominent feature of B — > Krj ) decays is the large 
disparity between the branching ratios for B — > Krj and 
B — > K-q decays. In the SCET framework this is quite 
naturally explained through a constructive and destruc- 
tive interference of different terms in the amplitudes as 
has been first suggested in |94l l95| . Specifically, the am- 
plitudes A B ^ K (') are related to Ab^k-q, and Ab^Kti s 
through a rotation (|41|l 



df) 



a b^k v > = cos Ms^Kns + sin Mb-Xi,, , ( 116 ) 
A B->R n = - sin Mb^k Vs + cos i>A s ^ R , (117) 



with cf> = (39.3 ± 1.0)°, so that cos<j> ~ sin<^. There 
is therefore a constructive interference in Ag^g-i and 
a destructive interference in A B _^ Rr) provided that 
Aq^kt] ~ A B^K-q s > which is exactly what is found 
in SCET. In B — > Krj ) there is a similar hierar- 
chy of terms that was found in B — > ttK decays so 
that charming penguin contributions are the largest, 



while tree contributions are A„ / A. 



0.04 sup- 



pressed, and QCD penguin and EWP contributions ~ 
max[C3, . . . , Cio]/as(2m c )Ci ~ 0.1 suppressed. The 



Ab-^k- =^rn%{x^ (A*** + A%°) +■■■} 

—Ago^ Vs K°J 



(118) 



AB°^rt q K°i 



A B -^ lqK - -^"^ r - c v ^ 



(119) 



where the ~ sign denotes equality up to smaller terms 
represented by ellipses. This gives for the Ab^k^ and 
Ab^Kti' amplitudes 

G ( 

A b -^ v k- =-^m B X { c s) cos0|(\/2 - tan0),4 cc , g 



G c 
A B -^r,'K- =—7=™^^ cosM (1 + v^tan^A 
V2 l- 



- tan (pj A cc H j , 



(120) 



cc,g 



1 + -!= t&nd) ) .1, 
V2 



}■ 



(121) 



where the SU(3) relations jS^I and (|54*|) were used. 
The same expressions hold for B° — > K°Tj( > ampli- 
tudes. The coefficient in front of A cc in the B — > Krj 
amplitude is (1 — v2 tan (j>)/{ v2 + tan (f>) = —0.07 sup- 
pressed compared to the one in the B — ► Krj' ampli- 
tude. The relative suppression of the A ccg contributions 
(V2 — tan0)/(l + V2tan0) = 0.28 is not as strong, so 
that B~ — > TjK~ and B° — > rjK are dominated by glu- 
onic charming penguin contributions. In the absence of 
the A ccg term, the corresponding branching ratios for 
B — ► rjK would be an order of magnitude smaller still, 
i.e. of order O(10" 7 ) instead of O(10~ 6 ). 

Note that the partial cancellation of charming penguin 
contributions in B — > rjK amplitudes occurs regardless 
of the strong phases carried by A cc and A cc g . In par- 
ticular, the strong phases of A cc and A cc ^ g that were ob- 
tained from the fit (H3, DHJl, 0HJ>, are near 180° and 
90° mod 180° respectively, so that there is little further 
cancellation between the two terms in <|120f) . A similar 
pattern of destructive and constructive interference has 
also been observed in the framework of QCD factoriza- 
tion mm. 

Note that the enhancement of BriB — > Krj') over 
BriB — ► ttK) is almost entirely due to the additional glu- 
onic charming penguin A ccg that arises from the annihila- 
tion of the n direction collinear quark with the spectator 
quark, where simultaneously two new n collinear gluons 
are emitted, Fig. These contributions correspond to 
singlet penguin amplitude s in the diagrammatic SU(3) 
approach, see appendix IO Explicitly, in the SU(3) limit 
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TABLE VIII: Predictions for the CP violating S parameters. The errors on the predictions are estimates of SU(3) breaking, 
l/m& corrections and errors due to SCET parameters, respectively. 



Mode 




Exp. 


Theory I 


Theory II 


B ^ 


ir r] 




-0.90 ±0.08 ±0.03 ±0.22 


-0.67 ±0.14 ±0.03 ±0.81 


B° — > 


ir° V ' 




-0.96 ±0.03 ±0.05 ±0.11 


-0.60 ±0.08 ±0.08 ±1.30 


B° — > 


IV 




-0.98 ±0.06 ±0.03 ±0.09 


-0.78 ±0.19 ±0.12 ±0.22 


B° — ► 


Wl' 




-0.82 ±0.02 ±0.04 ±0.77 


-0.71 ±0.14 ±0.19 ±0.29 


S ^ 


v'v' 




-0.59 ±0.05 ±0.08 ±1.10 


-0.78 ±0.09 ±0.19 ±0.23 


S ^ 
B° — > 


K s r,' 
K s n 


0.50 ±0.13 (S = 1.4) 


0.706 ± 0.005 ± 0.006 ± 0.003 
0.69 ±0.15 ±0.05 ±0.01 


0.715 ± 0.005 ± 0.008 ± 0.002 
0.79 ±0.14 ±0.04 ±0.01 



An- 



>ri'K- 



A* 



t+K- 



COS ( 



1.22+ 1.67- 



Ac 

V2 J A c 
\[2 sin 
A 



A 



ccy 



(122) 



ccy 



A,. 



where ellipses denote numerically smaller terms. In the 
limit Accg -> thus Br(B -> Krf) and Br(B -> irK) 
would be of similar size, while in SCET we expect A ccg ~ 
A cc which provides the observed enhancement. 

In QCD factorization, A ccg is perturbative and is pro- 



portional to the i<2 contribution in the F 



B — *rj,rj 



form 



factors of Ref. 13] . As previously discussed, Fo is not 
known from other sources and the authors of [13J as- 
sume the arbitrary values F2 — 0,0.1. Other mech- 
anisms that were proposed in the literature to explain 
the large Br(B — > Krf) are found to be either as{mb) 
or \jrrib suppressed in SCET. The contributions due to 
b — > sgg — > srj' coupling that would arise from integrating 
out the charm loop, Fig. 0,HEIi3 an d could be inter- 
preted as effective charm content of rf meson 0, |§H > 
lead to a 1/ ml suppressed operator (|14fl once matching to 
SCETn is performed. A mechanism in which one gluon 
is emitted from b or s quarks, with the other gluon com- 
ing from charm loop or from Og g insertion, Fig. |3 leads 
to either as (nib) or l/m& suppressed contributions as al- 
ready discussed below Eq. (fT^I . The har d sp ectator con- 
tribution b — > sg* g* — > srf discussed in |l00l |. where one 
of the hard off-shell gluons is emitted from the spectator 
quark, matches onto power suppressed SCETi operators 
with additional soft and collinear spectator quark fields 
obtained by integrating out the hard gluon (and other 
hard degrees of freedom). Similarly, the gluon conden- 
sate mechanism of Ref. |l0l| corresponds to a match- 
ing onto power suppressed operators with additional soft 
gluon fields once the intermediate hard-collinear gluon is 
integrated out. 

Since the smallness of Br(B — > r]K) arises from two 
large numbers cancelling, the predictions for this mode 
are rather uncertain, with modest variations on input 
parameters leading to larger relative variations on the 



observables. This could be used in the future to better 
constrain the SCET parameters C(J)gi A ccg . Of special 
interest are the direct CP asymmetries in B~ — > r\K~ 
and B° -> nK ° tha t can resolve between the two solu- 
tions (cf. Table EID • Defining the "tree" and ' penguin 
amplitudes as in 172fl . one has 



T B o_ 



Gf 

+ cf K 



2 hq 

(Ci 



C BK (c 2 
1 



1 

N 



N 



1 



(x-\„)Ci) 



T, 



(123) 



(124) 



where ellipses denote smaller terms coming from inser- 
tions of QCD penguin and EWP operators. The "tree" 
amplitude T B - ^ rjqK - receives contributions from con- 
figurations with two n-collinear gluons, Fig. that are 
part of SCET parameters (|43(l . In the diagrammatic 
approach these terms correspond to often neglected an- 
nihilation amplitudes as shown in appendix El We find 
them to be of LO in l/m& and should be kept in the 
analysis (in AS — 1 decays they are CKM suppressed 
and are thus numerically significant only for direct CP 
asymmetries). 

At LO in as(nib) the "tree" amplitudes Tgo^^go and 
Tb-^, u k- do not receive contributions from tree opera- 
tors 0" 2 > so that 

Tb ^j] s k<> ^ T B o^, h K<h T B -^ VeK - <C T B -^ r)qK - , 

(125) 

from which one obtains an approximate relation 



Tb^k T b ^ v , k 



cos ( 



sm i 



(126) 



where B (K) can be either B~(K~) or B° (K°). No 
such simple relation exists between P B ^ n ^ and P b ^„ik- 
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These are given by ((HO) and ((HTJ) (after division by A^ s) ) 
up to numerically smaller terms, and have a hierarchy 
P&^r)R ^ P&^ri'K as already discussed before. 
Defining 



-21m 



T, 



*W \ P, 



(127) 



B—*f 



where — 21-niXu^ / xi^ — 0.037, the CP asymmetries are 
to first order in small parameter rj 



Af = r f smS f + 0(rj). 



(128) 



The values of rf and the strong phase difference Sf be- 
tween "tree" and "penguin" amplitudes for B — ► rf'K 
decays are given in Table IIXI While the "tree" ampli- 
tudes in B — > rf ' K and B — > r\K are of approximately the 
same size (|126fl . the "penguin" amplitudes in B — > r/K 
are on the contrary suppressed as discussed above, lead- 
ing to an order of magnitude larger value for r g. Fur- 
thermore, in B — > rjK decays the strong phase differ- 
ence is predominantly determined by the gluonic charm- 
ing penguin A ccg due to a much larger suppression of 
A cc , so that the corresponding CP asymmetries are very 
sensitive to the value of arg A ccg . 

A different type of cancellation occurs in the B° — > 
Tr°r] amplitude. The two relevant amplitudes are (taking 



B — > Kr], the nongluonic contribution A™ is present 
only in Ago^^o^ and not in Ago^^a^ and therefore 

does not cancel in B° — > tt°t]. The same conclusions 
regarding charming penguins hold for B~ — ► n~r). Note 
that unlike B° — > Tr Q i] q , the "tree" term in A B -^ rjq ^ T - , on 
the other hand, is not predominantly gluonic since there 
is no equivalent cancellation to the one in the "tree" term 
of fTSPft . 

The predictions are fairly uncertain also for observables 
in B° — » rj^rf} decays, since these depend on both 

and C~, similarly to B — ► rf 'tt decays. We have 



Aga_> nv =A B0 ^ Vq11q cos cf> + A B(> ^ VaVa sin 

- A §o^ Ws sin 20, 
A- —(A- A- sSin2^ 

A B°^ijif —\ A B a ^ri q r lq — A B°^r ls vJ ^ 

+ 4b°-M7 9 »j. COS 20, 
A B°^ n ' v ' = A B«^ nq r h sm2 + 



+ A B°^ VqV , sin 20. 



(131) 
(132) 
(133) 



The amplitude Ago^,, „ receives contributions only 
from Qed,7d operators and thus has no charming pen- 
guin contributions. These are present in A B o^ Vqriq and 
■ A S°^n ris anc ^ therefore also in the amplitudes for the 
decays into mass eigenstates rf ^t/')- 



Am 



G 



B°^TT° 



F 2 r C 2 



„(<*) 



" s V2 " L V2 



+ •■•], 



(129) 



V2 



■^■B°^Tr°ri q 



G F 



V2 



2 

h (d) 



+ ^-(/.c;" 9 -/,,cr) 

-Ai rf) (A^ ff + ^)+... 



(130) 



where the ellipses represent the numerically smaller con- 
tributions from Q3d....,6d operators, while b% is evaluated 
as \l d) [C 2 + (1 + (x- 1 )* )Cx/N] + ■■■. Since f„ ~ U and 



+ 2C^ a , the tree amplitudes Tgo^o^') are 
predominantly coming from Cfpi" > i- e - f rom gluonic con- 



tributions. Furthermore, the combination c^Cs + 2 U; C/s 
is a linear combination of and C". Since the latter pa- 
rameter is largely unknown, cf. (|107fl . (|lllfl . this leads 
to relatively large errors on the predicted observables in 
B° — > 7T ?/ ) decays as can be seen from Tables IVTfl and 
IVIIII Measuring these observables would thus greatly 
improve our knowledge of £~ . 

The situation is reversed for charming penguin con- 
tributions. In B a — > 7r°?7 the contributions from A^ g 
partially cancel, just like in B — » _ftTr/. However, unlike 



C. S parameters in penguin dominated modes 

The CP violating S parameters in the AS = 1 de- 
cays B°{t) -> Ksxn , B°(t) -> Ks,lV (,) are of special 
interest because of the large CKM suppression of "tree" 
amplitudes over "penguin" amplitudes that was already 
discussed in the previous subsection. The decay ampli- 
tudes thus to a first approximation do not carry any weak 
phase and cancel in I|tj8fl . leading to the approximate re- 
lation Sf ~ ~ 7 lf F sm 2/3 where rff is the CP of the final 
state. This leads us to define an effective angle through 



Sf = 



-r^ p sin 2/3' 



/ 



(134) 



In this way the modes with K$ and Kl in the final state 
have the same (3 eS (neglecting CP violation in the kaon 
sector). 

In the decays at hand the common final state of B° 
and B° is provided by K° — K° mixing, that causes the 
mass eigenstates Ks,l to be an admixture of K° and K° 
states. For instance, we have 



A so ^ 



A 



B°^ V K S ,L 



Pk a b -+ v k° 



(135) 



with q K /pK = - V u S Vud/V us V* d ~ -1, where the phases 
of l^ ), \K°), \B°), \B°) were chosen so that CP\K°) = 
\K°) and CP\B°) = \B°). This relation can then be 
used in (|68|l to obtain the expression for Sf. As in the 
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TABLE IX: Predictions for "tree" over "penguin" ratios 77 
(first row of each mode, in terms of x 1CP 2 ) and strong phase 
differences 5f (second row), defined in Eq. 11271 . for several 
AS — 1 B° and B~ decays (separated by horizontal line). 
Theory I (II) column correspond to two sets of SCET param- 
eters Solution I (II) in itTUol - ltTTSll . Since B° -» tt°K° does 
not depend on isosinglet SCET parameters only one predic- 
tion is given. The meaning of errors is the same as in Table 

IVTTTI 



Mode 


Theory I 


Theory II 


rj'K 


2.9 ±0.7 ±0.7 ±0.4 


3.0 ±0.7 ±0.7 ±0.5 




(159 ±10 ±24 ±4)° 


(-117±11±24±6)° 


V K° 


21±20±4±3 


23 ± 24 ± 5 ± 4 




(101 ±58 ±20 ±4)° 


(-58 ± 62 ± 20 ± 7)° 


n°K° 


14±4±3±2 






(24 ± 20 ± 20 ± 6)° 




rj'K- 


2.8 ± 0.8 ± 0.5 ± 1.3 


0.8 ± 0.6 ± 0.1 ± 1.1 




(-22 ±11 ± 17 ±4)° 


(61 ± 11 ± 15 ±4)° 




34±31±7±3 


42 ± 46 ± 9 ± 5 




(105 ±57 ±20 ±3)° 


(-65 ±67 ±20 ±4)° 



previous subsections we decompose the amplitude into 
"tree" and "penguin" parts according to CKM element 

content Ag^j = Ac ± A« Tg_f . Due to the large 

CKM hierarchy in the AS = 1 decays "penguin" terms 
dominate over "tree", with the last term ~ A 2 = 0.04 
CKM suppressed compared to the first one, as discussed 
in Se ction IIII Bl Expanding in this small ratio we have 
1021 



AS f 



sin 2pf 



sin 2/3 



rf cos Sf cos 2/3 ± 0(rj), 



(136) 



where the "tree" over "penguin" ratio, r/e l5/ , was de- 
fined in Eq. I|127l) and already contains the ratio of CKM 
elements. 

The difference between sin2/3| ff and sin 2/3 vanishes 
if either Tf is zero or if the strong phase difference Sf 
between "tree" and "penguin" amplitudes is equal to 
±90°. The largest deviation, on the other hand, is ob- 
tained for Sf = 0°,180°. The deviation itself cannot 
be very large in the Standard Model because of the al- 
ready mentioned CKM suppression of order A 2 = 0.04, 
unless the ratio Tf/Pf is very large. This can happen 
in t\Ks,l decay modes, where there is large cancellation 
in -Pgo^jjA' l|120|) . No such cancellations are possible in 
rj ' K$,L modes. The large branching ratios of rj' K neces- 
sarily imply small deviations of sm2(3'ff Ks L from sin 2/3 
in the context of Standard Model. Using the values of 
parameters in l|86|) . (|87|l . and i|106|) - (|113[l obtained in the 



previous two sections 

As Th. { (-1.9 ± 0.5 ± 0.6 ± 0.3) x 10~ 2 
v ' Ks ' L ~ I (-1.0 ±0.5 ±0.8 ±0.2) x 10~ 2 

. „ Th. / (-3.4 ± 15.5 ± 



(137) 



__5.4±1.4)xl0- 2 , 
(7.0 ± 13.6 ±4.2 ± 1.1) x 10" 2 , 



where the upper (lower) rows correspond to Solution I 
(II) sets of SCET parameters in (fH)6>l(iT3|) . while 



AS 7r o Ks L = (7.7 ± 2.2 ± 1.8 ± 1.0) x 10" 



(139) 



where again the first two errors are due to SU(3) breaking 
and expected l/m^, a s (mb) corrections, while the last is 
due to experimental uncertainties on SCET parameters 
extracted from the x 2_n t to branching ratios and direct 
CP asymmetries. 

It is illuminating to evaluate the ratios r/e l<5/ Ijl27(l of 
"tree" and "penguin" terms for these modes. The nu- 
merical values are gathered in Table IIXI As expected 
from the general arguments outlined above, the ratio Tf 
is relatively large for B° — > t\Ks,l due to cancellations in 
"penguin" amplitudes. The corresponding strong phases 
Sf both for Solution I and Solution II sets of SCET pa- 
rameters (|106l) - (|113j) are not close to 0° or 180°, so the 
values of AS v k s l m P38|l are not close to the maxi- 
mal possible deviations for given values of 77. On the 
contrary, AS 7T o Ks L in l|139|) is already close to maximal 
positive deviation for fixed value of 77. Phenomcnolog- 
ically probably most interesting is ASti>k s l , whose ab- 
solute value is much smaller and is below 4% even if Sf 
is taken to be completely unknown. In accordance with 
general expectations very small values of | AS^Ks l I hare 
also been found in other ap proaches to two-body B de- 
cays; in QCD factorization |l4l 103|, in QCD factoriza- 
tion with modeled rescattering 104] and are con siste nt 
with bounds obtained using SU(3) symmetry [92lll05| . 

The predictions H138|l - l|139(l are to be contrasted with 
the experimental findings, where (negl ectin g the small 



difference between S 



AS*. 



ri'Ks 



Exp. 



and sin 2/3 ,106]) 



-0.23 ±0.13. 



and 



AS*. 



Exp. 



-0.41 ±0.26, 



(140) 



(141) 



while no information on AS^Ks IS Y e t available. The dif- 
ference between S^Ksi Stt k s an d sin 2/3 has been even 
more pronounced in the past, and has been reduced in 
the past year to the present level of almost 2 a. Since pre- 
dictions for these two quantities in SCET are not prone 
to large uncertainties as shown by the errors in (|137f) . 
(|139fl , a further reduction of experimental errors with un- 
changed central values would be a clear signal of beyond 
the Standard Model physics. 
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TABLE X: Predicted CP averaged branching ratios (xlO 6 , first row) and direct CP asymmetries (second row for each mode) 
for AS = and AS = 1 5° decays (separated by horizontal line). The columns Theory I (II) correspond to two sets of 
SCET parameters fTnSl - fTEfl . Since the decays into nonisosinglet mesons do not depend on parameters in 11061 - 11131 only one 
prediction is given. The errors on the predictions are estimates of SU(3) breaking, l/m;, corrections and errors due to SCET 
parameters, respectively. 



Mode 


Exp 


Theory I 


Theory II 


B° _» K-K+ 

5° 7V°K° 

50 r/0 
B]; —y rjK 

B s — > V A 


< 2.2/ d // a ° 

- 


4.9 ±1.2 ±1.3 ±0.3 
0.20 ±0.17 ±0.19 ±0.05 
0.76 ±0.26 ±0.27 ±0.17 
-0.58 ±0.39 ±0.39 ±0.13 
0.80 ± 0.48 ± 0.29 ± 0.18 
-0.56 ±0.46 ±0.14 ±0.06 
4.0 zt l.o zt U.4 zt U.o 
-0.14 ±0.07 ±0.16 ±0.02 


0.59 ± 0.34 ± 0.24 ± 0.15 
0.61 ±0.59 ±0.12 ±0.08 

o.y ± L.o ± U.o ± U.4 

0.37 ±0.08 ±0.14 ±0.04 


B° -> K~K+ 
5° -» K°K° 
B° -» W ° 


(9.5±2.0)/ d // s ° 


18.2 ± 6.7 ± 1.1 ±0.5 
-0.06 ±0.05 ±0.06 ±0.02 
17.7 ±6.6 ±0.5 ±0.6 
< 0.1 

0.014 ± 0.004 ± 0.005 ± 0.004 


0.016 ± 0.007 ± 0.005 ± 0.006 


B° - rf*» 




0.006 ± 0.003 ± 0.002to.oo6 


0.038 ± 0.013 ± 0.016±g:g|g 


fl° -» 7777 

s° - v'v' 




7.1 ±6.4 ±0.2 ±0.8 

0.079 ± 0.049 ± 0.027 ± 0.015 

24.0 ±13.6 ±1.4 ±2.7 

0.0004 ± 0.0014 ± 0.0039 ± 0.0043 

44.3 ± 19.7 ±2.3 ± 17.1 

0.009 ± 0.004 ± 0.006 ± 0.019 


6.4 ±6.3 ±0.1 ±0.7 
-0.011 ± 0.050 ± 0.039 ± 0.010 
23.8 ±13.2 ±1.6 ±2.9 
0.023 ± 0.009 ± 0.008 ± 0.076 
49.4 ± 20.6 ± 8.4 ± 16.2 
-0.037 ± 0.010 ± 0.012 ± 0.056 



"The production fraction ratio of B 1 ^ mesons is / d / f s 4 60] . 



D. B s decays 

Using SU(3) symmetry allows us to make predictions 
for B® decays as well. Predictions made with the values 
of SCET parameters ®, and l(lH5 )l -l(TT5 )l for CP 

averaged branching ratios and direct CP asymmetries are 
collected in Table Ixl while the predictions for the observ- 
ables (Sf) B , © and {H f ) B , fl® are given in Table IXll 
The SU(3) breaking on the SCET parameters relations 
(HHJ-HEJ) and was assumed to be 20% with a 20° 
variation on the charming penguin's strong phases. The 
second errors in Tables |X] IXII estimate the remaining 
order 1/rrtf, and as{mi) corrections. These are obtained 
from a 20% variation on the size and a 20° variation on 
the strong phase of the leading order amplitudes propor- 
tional to Al /) or \[ f) . 

Many observations made about B° and B~ decays 
hold also for B° decays. For instance AS = 1 decays 
B® — ► KK and B® — ► rf'rf ' are dominated by nonper- 
turbative charming penguins due to a CKM hierarchy 
just like B — > nK, B — * Krf ^ decays. Expanding in the 
CKM suppressed "tree" over "penguin" ratio rt <|127|) 



the observables from time dependent decays (|69(l . (|70ll 

(S f ) Ba = ?// P sin 2e-r)f p r f cos6 f cos2e + 0(rj), (142) 
and 



(Hf)i 



vf P cos 2< ! 1 - 



t?j p sin 2e 



Re 



Im 



(AgAg) 



( cos 2 5f — 



Tf COS Sf 

(143) 



while the expression for direct CP asymmetry to first 
order in rt is given in () 128(1 . Since e ~ 1° in the Stan- 
dard Model, (Hj)b s for the penguin dominated decays 
B® — > KK and —* rf'rf} is expected to be very close 
to 1. In the Standard Model sin2e ~ 0.035 and thus 
sin2e ~ r/, so that the deviations of [Hf) Bs from unity 
are numerically of order 0(r 2 *). For 1 ^> 77 > 2sin2e, 
the r 2 correction in the first term in (|143fl is actually 
larger than the second term in 1)1431) which starts at lin- 
ear order in rj. In this case rfi \Hf) Be is smaller than 
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TABLE XI: Predictions for (S'/)s 3 (first row in each mode) and {Hf)B s (second row) parameters in B s decays. The columns 
Theory I (II) correspond to two sets of SCET parameters (f lOfipi - (|1 lHp> . Since the decays into nonisosinglet mesons do not depend 
on parameters in 11061 - 11131 1 only one prediction is given. The errors on the predictions are estimates of SU(3) breaking, 1/mb 
corrections and errors due to SCET parameters, respectively. No predictions are made for (S/)b s and {Hf)a s in B° — > K°K° 
and B® — > it°r)' , see text. 



Mode 


Exp 


Theory I 


Theory II 


B° a -» K S n° 


- 


-0.16 ±0.41 ±0.33 ±0.17 






- 


0.80 ±0.27 ±0.25 ±0.11 




B° -> Ksv 


- 


0.82 ±0.32 ±0.11 ±0.04 


0.63 ±0.61 ±0.16 ±0.08 




- 


0.07 ±0.56 ±0.17 ±0.05 


0.49 ±0.68 ±0.21 ±0.03 


B° - Kgrf 


- 


0.38 ±0.08 ±0.10 ±0.04 


0.24 ±0.09 ±0.15 ±0.05 




- 


-0.92 ±0.04 ±0.04 ±0.02 


-0.90 ±0.05 ±0.05 ±0.03 


B° -> K~K+ 


- 


0.19 ±0.04 ±0.04 ±0.01 








1 - (0.021 ± 0.008 ± 0.007 ± 0.002) 




B° a -> A 




0.45 ±0.14 ±0.42 ±0.30 


0.38 ±0.20 ±0.42 ±0.37 






-0.89 ±0.07 ±0.21 ±0.15 


-0.92 ±0.08 ±0.17 ±0.15 


b°s - m 




-0.026 ± 0.040 ± 0.030 ± 0.014 


-0.077 ± 0.061 ± 0.022 ± 0.026 






1 - (0.0035 ± 0.0041 ± 0.0019 ± 0.0015) 


1 - (0.0030 ± 0.0048 ± 0.0017 ± 0.0021) 






0.041 ± 0.004 ± 0.002 ± 0.051 


0.015 ± 0.010 ± 0.008 ± 0.069 






1 - (0.0008 ± 0.0002 ± 0.0001 ± 0.0021) 


1 - (0.0004 ± 0.0003 ± 0.0003 ± 0.0007) 


B° - VV' 




0.049 ± 0.005 ± 0.005 ± 0.031 


0.051 ± 0.009 ± 0.017 ± 0.039 






1 - (0.0012 ± 0.0003 ± 0.0002 ± 0.0017) 


1 - (0.0020 ± 0.0007 ± 0.0009 ± 0.0041) 



1 irrespective of the strong phase difference 5f. The de- 
viations of {Hf)s s from 1 for different modes are listed 
in Table EJ 

The charming penguin dominance in the KK modes 



A m _ 



Gj 



>K-K+ 



V2 

-Am^KORo 



KK 



(144) 



furthermore leads to an approximate relation 

r s°— k-k+ — ^b»^k°k°- (145) 

The amplitude for B® — > K°K° does not receive any LO 
contributions from Qi] 2 s operators, cf. Table IIVI so 
that at this order in the as (tub) expansion the "tree" 
amplitude is entirely due to QCD penguin and EWP 
operator insertions. At NLO in as(mg) there could, 
however, be contributions proportional to Ci^otsijnB) 
leading to a large correction to the predicted direct CP 
asymmetry. We thus give in Tabled only an estimated 
upper bound of this observable. Similarly, no prediction 
of (Sf)B, , (Hf) b s for this mode is made in Table IXll but 
conservatively we can expect 



In B2 



\(Sk s k s )b 3 \ < \(S k + k -)b 3 I, 
(Hk s Ks)b 3 \ < |1 — {Hk+k-)b 3 



(146) 
(147) 



r, y "rf > amplitudes the dominating contribu- 



tions are also given by charming penguins 
G F 



Vs Vq 



Ab°-H7„tj. 



V2 B 
Gf 2 
V2 



(148) 



AW2(4«-( S )+4«( S )) + 



(149) 



with the ellipsis denoting smaller contributions, while in 
Ag°^ri j) there are no charming penguins. Using l|131|) - 
(|133[) we then get for the coefficients in front of A cc , A ccg 
in the SU(3) limit (apart from the common multiplicative 
factor GfTOb/v 7 ^) 



mode 


A C c 


A 

-^ccg 




-> 7777 


\/2sin 2 


V2sin 2 


<f> — sin 2cf) 


m 


— > 7777' 


sin 20 


cos 2(f) — 


sin 20 


m 


— > 77' 77' 


a/2 cos 2 4> 


\/2cos 2 


<f> ± sin 2(f) 



which for <p = 39° gives numerically 



mode 


Acc 


Accg 


m 


-> 7777 


-0.42 


0.56 




— > 7777' 


-0.48 


-0.69 


m 




1.83 


0.85 



explaining the pattern of branching ratios in Table El 
Note that the branching ratios for B® — * rjrj and B® — > 
77' 77' have an additional symmetry factor 1/2 in (|65|l be- 
cause of indistinguishable particles in the final state. 
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TABLE XII: Predictions for "tree" over "penguin" ratios 77 
(first row of each mode, in terms of x 10~ 2 ) and strong phase 
differences 5/ (second row), defined in Eq. itHTl . for pen- 
guin dominated AS = 1 B° decays. Theory I (II) column 
correspond to two sets of SCET parameters Solution I (II) in 
11U6I I- 1113II . Since B° — > K + K~ does not depend on isos- 
inglet SCET parameters only one prediction is given. The 
meaning of errors is the same as in Table IVTrTl 



Mode 



Theory I 



Theory II 



K~K+ 
'PI 



rfrf 



16.2 ± 3.5 ± 3.3 ± 1.0 
(-159 ±18 ±21 ±5)° 
10 ± 5± 2± 2 
(53 ±25 ±20 ±4)° 
0.5 ±0.4 ±0.2 ±5.1 
(175 ± 17 ±48 ±4)° 
1.6 ±0.5 ±0.4 ±3.6 
(145 ± 13 ±25 ±9)° 



12±7±2±2 
(-6 ±26 ±20 ±4)° 
3.1 ± 1.0 ±0.5 ±0.5 
(47 ± 18 ± 18 ±4)° 
4.0 ± 1.1 ± 1.0 ±3.6 
(—112 ± 12 ±25 ±9)° 



Since B s — > KK and B s — ► rf^rf^ decays are penguin 
dominated, they represent ideal probes for new physics 
searches. A very useful observable in this respect is 
(Hf)s, (|143fl . I n the Standard Model we have a very 
robust prediction that (-ff/)s 3 — 1 for penguin domi- 
nated modes B s — > KK and B s — > 7/ '7/ \ where the 
deviations from this relation are of order O(r^) as dis- 
cussed above, where 77 itself is of order A 2 ~ 0.04 (|127fl . 
The predicted values in SCET for different modes are 
given in Table IXlTl The deviations from (ii/)s 3 = 1 are 
therefore expected to be very small, well below percent 
level for B s 



jy( )^( ) modes. If virtual corrections from 
the beyond Standard Model particles modify either the 
B s — B s mixing phase or introduce new phases in b — > s 
penguins, a hint of which may have already been experi- 
mentally seen in B — > ttK , B — > rf K modes, (-ff/)s 3 can 
easily deviate from 1 by a correction of 0(1). Another 
attractive feature of ( Hf) r „ is that it can be measured 
in untagged B s decays iti-'j . 

A special case among the AS" = 1 B s decays are 
the decays B® — ► 7r ?7^ '. Since these are AI = 1 
transitions, there are no QCD penguin contributions to 
these decays, including no contributions from the non- 
perturbative charming penguins. The theoretical con- 
trol over these modes is thus much greater, on par with 
B~ — » 7T~7r°, where charming penguins are also absent. 
Furthermore, B® — + tt°t]( ) decays are color suppressed, 
which makes them a perfect probe of the SCET predic- 
tion that color suppression is lifted due to Q^J contri- 
butions 
order ~ 
of the " 
in B° s - 
(Sfh. 



Unfortunately their branching ratios are of the 
*> 10~ 7 — 10~ 8 because of the CKM suppression 
"tree" amplitude. The absence of QCD penguins 
— * TT ^ ) also leads to a generic prediction that 
and 1 — (Hj)b s can be of 0(1). Namely, the 
"penguin" amplitude in B° s — ► 7r°77^ ' is coming exclu- 
sively from EWP operators, which have small Wilson co- 
efficients, making them of similar size as the CKM sup- 



pressed "tree" amplitudes, giving rt ~ 0(1) in (|142fl and 
(|143[) . Furthermore, because there are no contributions 
from charming penguins, while the other contributions 
factorize at LO in the strong phase differences 

vanish at leading order in a s (rrib) and l/m,b, 5f = (0° 
mod 180°), so that A^ p ((Hgj) in B° s -> ttV'' vanishes 
at this order. This will be lifted at NLO. 

Scanning the allowed range of the at present very 
poorly constrained parameter £~ i|107|) . i|lll|) . we find 
that a cancellation between different terms in Ago^^a^ 
is possible for special values of C~. We thus give only 
a la range for the decay width, while no prediction on 
(Sf)B s and 1 — {Hf)B s for B® — > ir°r}' can be made. 

The AS" = decays have "tree" and "penguin" contri- 
butions of similar size. Furthermore, both B® — > K ^ ) 
and B® — > ir~K + , ir°K° do receive contributions from 
nonperturbative charming penguins, leading to sizeable 
direct CP asymmetries. A peculiar case is the B® — » 
K°r] decay, where a similar cancellation occurs as in 
B° — > K°rj between charming penguin contributions due 
to 77— 77' mixing. Namely, up to numerically smaller terms 
from Qzd,...M insertions we have in the SU(3) limit 



a _ _ Gf 2 



G 
72 



F 2 



\^{A CC + A ccg ) 



\^(-j=A cc + V2A c 



j_ hi (fid) , r ,(d)\ 

\/2 



(150) 



(151) 



giving a (1- v / 2tan</.)/(\/2±tan0) = -0.07 suppressed 
coefficient in front of A cc in B® — » K°rj compared to 
the coefficient in front of A cc in B® — > K°r]' . Simi- 
larly, the coefficient in front of A ccg in B° s — > K°r) is 
(1 — tan0/v / 2)/(l/v / 2 ± tan0) = 0.28 suppressed com- 
pared to the one in B® — > K°rj' . No such cancellations 
are present in "tree" amplitudes. Because B® — > K ^ ) 
decays are AS = 0, with no CKM hierarchy between 
"tree" and "penguin" amplitudes, the effect of cancella- 
tions on predicted branching ratios is less drastic than in 
B° — > K°r)( ), as can be seen from results in Table Ixl 



IV. CONCLUSIONS 

We have provided expressions for decay amplitudes of 
B°, B~ and B® mesons to two light pseudoscalar or vec- 
tor mesons, including isosinglet mesons r},r)',u},4>, using 
Soft Collinear Effective Theory at LO in l/m&. For the 
decays into 77, 77' mesons the contributions where a spec- 
tator quark is annihilated by the weak operator, while 
two collinear gluons are created, is found to be of the 
leading order in l/rn&. This leads to a new jet function 
in SCETi — > SCETn matching for the operators that 
contribute only to isosinglet final states. 

In the phenomenological analysis we work at LO 
in as (77ib) and 1/mf,. In particular, no expansion in 
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as(v / Amb) at intermediate scale is made. Instead, fol- 
lowing Refs. 0, El a set of new nonperturbative SCET 
functions is introduced, which are then extracted from 
data using a x 2 -fit. The numerical analysis is performed 
for the decays into two pseudoscalars. From separate x 2 - 
fits to B — > tttt,ttK data and to B — ► ■n-rf- \Krf ■* data 
SCET parameters for nonisosinglet final states and for 
isosinglet final states are extracted, respectively. Due to 
scarce data for isosinglet final states, SU(3) symmetry 
is imposed on the SCET parameters. Even so, two sets 
of solutions are found for the SCET parameters C(j)gi 
A ccg that enter only in the amplitudes for decays into 
isosinglet final states. 

Predictions for branching ratios, direct and indirect CP 
asymmetries in B°, B~ and B® decays together with es- 
timates of theoretical errors due to SU(3) breaking, 1/mj 
and as{rrib) corrections and errors due to extracted val- 
ues of SCET parameters are listed in Tables IVIVIIII and 
Tables IXlXll We find: 

• A discrepancy between theoretical and experimen- 
tal values for the R and R n ratios of decay widths 
(JHU, HUl, EDI), §12} at a level of 2a if both tttt 
and irK data are used simultaneously to determine 
the SCET parameters. No statistically significant 
deviations are found if only a fit to ttK data is 
made, but the resulting values of SCET parame- 
ters strongly violate flavor SU(3). 

• The enhancement of Br(B — > Kr/') over Br(_B — > 
Krf) is naturally explained in SCET due to destruc- 
tive and constructive interference governed by rj—rf 
mixing as first proposed by Lipkin [94, 95] . The en- 
hancement of Br(_B — > K-q') over Br(£? — * irK) is 
due to the "gluonic" charming penguin A ccg , where 
the n collinear quark annihilates with the spectator 
quark producing two n collinear gluons. The other 
mechanisms for enhancing Br(_B — ► Krf) discussed 
in the literature are found to be suppressed by at 
least one order of l/mb- 

• The SCET parameters C(J)g an( i A ccg correspond- 
ing to annihilation diagrams with two simultane- 
ously emitted collinear gluons are found to be both 
parametrically and numerically of the same or- 
der as the other leading order contributions in ac- 
cordance with SCET counting. The combination 
Qj g — £ g is at present very poorly constrained, but 
can be constrained by future measurements of de- 
cay widths and CP asymmetries in B — > 7r ?/ - 1 and 
CP asymmetries in B — > K°r]( > decays. The or- 
thogonal combination Q g + ( g which also enters 
into B — > rf > form factors is already much bet- 
ter contraincd and the corresponding values of the 
B — ► ri q _ s form factors are given in Eqs. (|114|l and 
QI 151) . In the future, with more abundant data on 
isosinglet decays, the assumption of SU(3) flavor 
symmetry that was used in the present analysis can 
be relaxed. This will significantly reduce the the- 



oretical uncertainties associated with many of our 
predictions. 

• The deviation AS of the S parameter from sin 2/3 
in B°(t) — > Ksrj' decay, which was not used as a 
constraint in the x 2_n t? is predicted to lie in the 
range [—0.026, 0] at la and to be below 4% even if 
no information on the strong phases is used. Fur- 
ther constraints on AS in B° — ► Kgr], Ksir decays 
are given in subsection IIII CI 

• Predictions are made for the branching ratios in 
B s — ► PP decays, as well as for related direct 
CP asymmetries, (S/)b s parameters and the co- 
efficient (Hf)s s multiplying sinh(Ari/r) in the 
time dependent decay width. A robust prediction, 
{Hj)b, — 1, for penguin dominated AS = 1 de- 
cays holds in the Standard Model up to corrections 
that are at the permil level for B s — > rfi >rf ** de- 
cays and at the percent level for B s — > KK decays. 
This prediction follows from general arguments in- 
dependent of the SCET framework. In beyond the 
Standard Model scenarios, the relation (-ff/)s s = 1 
can generically receive corrections of 0(1), making 
them very interesting probes of new physics. 
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APPENDIX A: NOTATION 

Choosing n» = (1,0,0,-1) and tt/ 1 = (1,0,0,1), so 
that fi = 7° + 7 3 and ^ = (7 — 7 3 ), while in shorthand 
notation n-k = k + and fi-k = k— for any four- vector fc^, 
and choosing for the Levi-Civita tensor normalization to 
be 

e 0123 = -1, (Al) 

while 75 = Z7°7 1 7 2 7 3 , so that for instance 
Tr[757 M 7 I/ 7 a 7' 3 ] = 4ie MlyQ ' 3 , and furthermore defin- 
ing 

ef = \e^n a n p (A2) 

so that e x l = -e^ 1 = 1 and 75 = §^,?i] 6^7777^, this 
leads to the following relations 

rttf±.PL,R = (9T ± i*T ) *Pl,r, (A3) 
WtfiPL.R = (sr T ief ) m,R, (A4) 

where g^ = g^ ' -\(n^n h ' +n^n v ) and P L>R = (lT7 5 )/2. 
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APPENDIX B: THE STRUCTURE OF JET 
FUNCTIONS 

In this appendix we show that in the matching l|24(l 
only four jet functions appear to all orders in as{yAmb)- 
Let us first discuss the case where in 

T[(t n WUigBl% u P R ,il ia (0)[igPiwH n ] jb (y), (Bl) 

the gluon fields are contracted, leading to J± and J jet 
functions l|24|) . At higher orders in as(v / A™f>) the cor- 
rections will come from insertions of leading order SCET 
Lagrangian 

4? = Zn, P > [in -D c + z PcT^-J P^j \tn, v , (B2) 

and the purely gluonic Leg ■ Thus the most general form 
of the operator with two external collinear quark fields 
that 1|B1[) matches onto is (up to operators with nontrivial 
color structure) 

S ab [< 1[(^ ± 'jT +1 i)P R ,L\ 1 [li X] 3 , 

i 3 

(B3) 

with 7^ coming from jS^, ft coming from the collinear 
quark propagator, and pairs of 7_l matrices coming from 
C% . After Fierz transformation then (up to an overall 
factor) 

1 3 (B4) 

i 

where 

T/ <8> T'j = ft® ft- ^7 5 <K> ^75- ffcf±® flj±u = , . 
= ^(1=F7 5 )0^-^7±<H)^7-L ! /. 

In the last equality the action of 75 on Pn,L has been 
used. The two Dirac structures give rise to the operators 
with jet functions J and J± in Q24JI. For instance the 
=F 75)® yl Dirac structure leads to the operator 

[q n , X u,m,R(/ n< -^] (fijlP R ,L,Y\ (B6) 

which follows from the relation fi fi fi oc and the fact 
that all the Lorentz indices in l|B4|) are contracted pair- 
wise (contraction with external p± momentum is sub- 
leading) except for the index a which is carried by the 
remaining 71. Namely, the 7j_ matrices with contracted 
indices can be permuted to be next to each other so that 
the end result is a pure number times 7". 

The /l7^<8> /i7_Li/ Dirac structure, on the other hand, 
leads to the operator 

[qn, IU jfrydVn, -xj (^Lllllf- (B7) 



This immediately follows from the relation [odd] 7^ [even] 
oc 7^7^ to be proven below. Here [odd] and [even] denote 
products of odd and even number of 7_l matrices, with 
7" either in [odd] or [even], while all other indices are 
contracted. 

Before we proceed let us show by induction that 
[odd]7^[odd]=0, where the indices of 71 matrices in 
[odd] are all contracted pairwise. This is true at low- 
est order, 7^7^7j_/3 = 0. Now assume that the rela- 
tion holds for A — 2 71 matrices and look at the case 
of A 7j_ matrices. Since in [odd] there is an odd num- 
ber of 7j_ there must be at least one pair of 7j_ matri- 
ces that has contracted indices and sits on the opposite 
sides of 7j\ Moving these matrices next to 7^ using 
7^7^ = 2^-7j_7j leads to terms of form [odd] 7 ^[odd] 
with A — 2 matrices, which are zero by assumption, and a 
term [even] 7^ 7^ 71/3 [even], which is also zero. Similarly, 
one can show that [even] 7^ [even] oc 7^. This holds triv- 
ially at lowest order with [even] empty. Let us assume 
that it holds for A — 2 j± matrices and move to A 7_l ma- 
trices. If in [even] 7^ [even] there are no cross contractions 
then [even] is just a number. If there are cross contrac- 
tions, then as before the corresponding two matrices can 
be moved next to 7^. This leads to terms [even] 7^ [even] 
with A — 2 matrices, which are proportional to 7^ by as- 
sumption, and a term [odd] 7^7^7.1/3 [odd], which is zero. 

We can now show by induction that [odd] 7^ [even] oc 
7^7^. The relation is trivially satisfied when [even] is an 
empty set and [odd]= 7" . Let us assume that the relation 
holds also for A — 2 7j_ in [odd] and [even] and move to A 
7j_. We distinguish two cases, (i) 7" is in [odd] and (ii) 
7" is in [even] . For option (i) the matrix 7" can be moved 
to the far left using 7^*7" = 2<?^ tQ — 7"7^ 1 - The terms 
with g^ a have A — 2 71 and are proportional to 7^7^ 
by assumption, leaving a term 7" [even] 7^ [even] oc 7^7^ 
(since [even] 7^ [even] oc 7^). In the case (ii) is moved 
to the right. This leads to terms with A — 2 gamma 
matrices, which are proportional to 7^7^ by assumption, 
and a term [odd] 7^ [odd]7j , which is zero as shown above. 

Finally let us discuss the matching in I|B1|) with 
fermion fields contracted. Generally, this leads to an op- 
erator of the form 

S ab bl i U^±< i+1 ^P^fnt^li-^ (B8) 

i 

where one of the Lorentz indices is equal to a, while the 
others are contracted pairwise. This general operator can 
always be written as a linear combination of operators 

S ab hi 7lP R ,L} ji Bt- Su ■ Bt_ xu , (B9) 
S ab [Bt_ £w ^P R , L fB^_ xu , (BIO) 
S ab [Pin- x ^PR,L] ji Bt«_ Sui . (Bll) 

This is trivially satisfied for only 7^ in (|B8(1 without ad- 
ditional 7j_ insertions. It is also true in general, since 
the Dirac structure [any] 7^ [any] 7^ [any] 7" [any], where 
[any] is a product of an arbitrary number of 71 matrices, 
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with the unshown Lorentz indices contracted, is a sum 
of 7^7^7" and the similar terms with permuted Lorentz 
indices. Decomposition in terms of operators ljB9|) - (|Bll|) 
then follows from 



v u ql li>v a 

7_l7±7_l = 91 7± 



(B12) 



Since (|B10(I and (|B 1 1|) differ merely by x <-> x inter- 
change, only one of the two is needed after integration 
over x in (|24|l . 



APPENDIX C: SU(3) DECOMPOSITION 

In the limit m s , Tn u .d <CAa useful approach is to use 
the transformation properties of the weak Hamiltonian 
@ under flavor SU(3) and decompose the d ecay ampli- 
tudes in terms of reduced matri x element s [107]. The 
diagrammatic approach of Refs. |l08l Il09j is equivalent 
to the SU(3) decompositio n of t he am plitudes (for recent 
applications see e.g. [30l IllCt llllj l. It is usually fol- 
lowed, however, by further dynamical assumpt ions, with 
annihilation and exchange topologies neglected |llOj . For 
nonisosinglet final states this assumption can be justified 
using SCET, since the reduced matrix elements corre- 
spondingto the two topologies are found to be 1 jm^ sup- 
pressed H 113, EH . As we will show in this appendix, 
the annihilation and exchange topologies, on the other 
hand, cannot be neglected for isosinglet final states since 
they come as leading order contributions in l/mf, ex- 
pansion in SCET. In this appendix we also provide the 
translation between our LO SCET results (|37J) and the 
diagrammatic language. All the results will be given as- 
suming exact SU(3), using thus the relations (|45|l . (|46|l . 
(|53(l , Il54|) along with a similar relation for the decay con- 
stants f M = U = fl< = f Vq = fr h ■ 

The effective weak Hamiltonian (|2l transforms under 
flavor SU(3) as 3®3<g)3 = 3©3©6©15, so that it can be 
decomposed in terms of a vector if* (3), a traceless tensor 
antisymmetric in upper indices, if]/ 7 '' (6), and a traceless 
tensor symmetric in the upper indices, ff] 1 ^ (15). We 
further define a vector of B fields Bi = (B + , B°, B®) and 
a matrix of light pseudoscalar fields 



/si _ M. 

/ J2 Jfi 



V2 Ve 

7T+ 



n 

.si _ M. 
K° 



K° 



+ ^t, ( ci) 



where the SU(3) singlet T70 ~ (uu + dd + ss) / v3 and 
octet r/s ~ (2ss—uu—dd)/\/6 admixtures of f? 9iS are used, 
which is the natural choice in unbroken SU(3). The most 
general Hamiltonian that has the same transformation 
properties under SU(3) as the weak Hamiltonian in @ 
is then a sum of terms that contribute to both SU(3) 



singlet and octet final states 

A^B i H i (2>)M 3 k M k + C 3 B % MiMlH k {3) 
+A 6 B l H l k J (6)M l J Mf + C 6 BiM)Hi k {%)M l k 



(C2) 



A 15 B t Hl 3 (15)MjMf + C lb BiM]Hf K (15) M, 



k ■ 



and terms that are nonzero only if SU(3) singlet is in the 
final state 



y : .H,M[iri:UMl + V 3 B l H l (3)M J j M k 
+V 6 B l Hl j (6)M^M l l + D 15 BiH% (15)M*M/, 

where for AS = decays 



(C3) 



if 2 (3) = 1, (C4) 
ffi 2 (6) = -fff(6) = if 3 23 (6) = -fff(6) = 1, (C5) 



2ifi 2 (15) = 2if 21 (I5) = -3fi 22 (15) = 

= -6if 23 (15) = -6fff(15) = 6, 



(C6) 



with the remaining entries zero, while for AS = 1 de- 
cays the nonzero entries in ff*(3), H k (&), H l k J (15) are 
obtained from (|C4|) - (|C6() with the replacement 2^3. 
Note that since the final state \PP) is symmetric, ther e 
are only 9 reduced matrix elements in B — ► PP |107| . 
Namely, the coefficients Cq,Aq,T>q in the above decom- 
positions (|C2|) . I|C3(1 always appear in the combinations 
Ce — Aq and T> e + A§. 

In the diagrammatic approach the linear combinations 
of reduced matrix elements £3,6,15 are redefined as t, c,p 
amplitudes 



< = 2C 6 + 4Ci 5 , 
c= -2C 6 + 4Ci 5 , 
P = C 3 — C(, — C15, 



(C7) 
(C8) 
(C9) 



while ^3,6,15j leading to amplitudes e,a,pa in the di- 
agrammatic notation, arc usually neglected. This dy- 
namical assumption can be justified usi ng SCET, where 
-^3,6,15 are found to be 1/mj, suppressed |9 Ulfllil"2| . while 
the remaining amplitudes are at LO in l/m& and as(rrib) 



G ^™2 , \ (h (d) h (dhf,,(d)_(dh 



(CIO) 



c = ^m|/ M [(4 rfJ - TO + (4 dJ - 4 d) )C], (Cll) 



p = ^=m% [f M (bi d) tj + eft) + \^A CC ] . (C12) 

The Wilson coefficients bf are here understood to be 
already convoluted with light pseudoscalar LCDA. In 
the SU(3) limit this amounts to a replacement mb/t02 = 
-m b /u! 3 = (x^ 1 )^ = {x^) K = (x^ 1 ) v ~ 3 in Eq. |JT§J. 
The amplitudes for AS — 1 transitions are obtained from 

^m-RM with -> b^, c^. 

The complete SU(3) decomposition of amp litud es for 
decays not containing 770 is given in Ref. [l09j both 
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TABLE XIII: The SU(3) decomposition of AS = (above horizontal line) and AS = 1 decays (below horizontal line) into final 
states with rjo- Each amplitude should be divided by the common denominator in the "Factor" column, so that for instance 
a b-^- vo = ( 2C 3 + C 6 +3Ci 5 + 2Ae+ 6Ais + 3£ 3 + 3£> 6 + 9£>i 5 )/\/3 = {c + t + 2p + t s +p s )/V3. The "Diagrammatic" column 
shows the decomposition in the diagrammatic approach notation, with .43,6,15 neglected. 



Mode 


C3 


C6 


Cis 


1 


1 

Ae 




£3 


£>3 


T>6 


T>15 


Factor 


Diagrammatic 


B~ — ► tv^t/o 


2 


1 


3 





2 


6 


3 





3 


9 


V3 


f. + c + 2p + t s + p s 


fso 
B" -> 7T u r? 


-2 


— 1 


5 





-2 


10 


-3 





-3 


15 


V6 


—2p + c s —p s 


B° -> r/ s rj 


-2 


3 


-3 





6 


-6 


-3 





9 


-9 


3\/2 


-(2c + 2p + c s +p s ) 


B° -> r/ rio 


2 








6 








6 


18 








3 


2(c + p + c s +p 3 + s ) 


B° s -» K° m 


2 


-1 


-1 





-2 


-2 


3 





-3 


-3 


a/3 


c+2p + p s 


B~ -» K-T/o 


2 


1 


3 





2 


6 


3 





3 


9 


a/3 


t' + c + 2p' + t' s + p' s 


B° -» tf% 


2 


-1 


-1 





-2 


-2 


3 





-3 


-3 


a/3 


c' + 2p'+p' s 


B s ° - Ao 





-2 


4 





-4 


8 








-6 


12 


a/6 


c' + c' s 


S° -> ??8?70 


4 





-6 








-12 


6 








-18 


3\/2 


c > + 4 p > - c > g+ 2p ' s 


5° -> 770 ??0 


2 








6 








6 


18 








3 


2(c' + p' + c' s + p' s + s' ) 



in terms of reduced matrix elements -43,6,15 > £3,6,15 as 
well as in terms of diagrammatic amplit udes and will 
thus not be repeated here. (In Ref. |l09j a differ- 
ent phase conventions was used, so the replacements 
7T° — ► — 7r°,7r~ — > —ir~,K~ — » — K~ need to be made, 
while the amplitudes into two indistinguishable states 
should be multiplied by \/2 to have our normalization 
of the amplitudes. In addition we find the "Factor" in 
Table 2. of [To3 for B s n r] 8 to be -\/3.) The complete 
SU(3) decomposition of amplitudes for decays into SU(3) 
singlets, on the other hand, is provided in Table IXlTTl 

The reduced matrix elements £3 and 2?3,6,i5 describ- 
ing the decays into SU(3) singlet final states are found 
in SCET to be all nonzero already at leading order in 
l/mft and as (mi,). Defining the singlet "diagrammatic" 
amplitudes 

t s = 6(P 6 + 2P 15 ), (C13) 
c s = -6(V 6 - 2D 15 ), (C14) 
Vs = 3(C 6 -C 15 -V 6 -V 15 +£ 3 ), (C15) 
s = 9(V 3 +V (i -V 15 ), (C16) 

we find at leading order in 1/mb and as(jrib) 
t s =^m%3f M [b^O g +cf( g ], 
c s =^m%3f M [(bf - b^)O g + (4 d) - 4 d) )C s ], 



p s =^m%3[f M ((bf - bf )(j + (cf - 4 d) )C) 



+ fM{b { f ) (j g + c[ d) ( g )+\WA ccg ] : 
Gp 

7T 



s o — ; 7^ m B^JM L(o s -o 6 )C/9 + (c 5 -c 6 )C, 



(C17) 



where so contributes only to 770^70 decays. The AS" = 1 
amplitudes t' s ,c' s ,p' s , s' a are obtained by replacing d — > s 
in l!C17|l . Note that these amplitudes are arising from 
gluon content of isosinglet final states and correspond to 
diagrams on Fig. 01 b) and0]d) and on Fig. with p s 
receiving also nongluonic contributions. 

In the applications of SU(3) decomposition using the 
diagrammatic approach it is frequently assumed that 
only one reduced matrix element, £3, is nonzero , while 
f 3,6,15 are assumed to be suppressed |lf)9l IllOj ]. This 
corresponds to taking p s nonzero, while neglecting t s , c s 
and so (commonly s = p s /3 is introduced instead of p s ). 
In the SCET result (|C17I) this would correspond to a 
limit C(J)g ^ C(J)> while from SCET scaling of the oper- 
ators we expect Oj)g ~ C(J> A-t present both possibilities 
are still allowed by experiment, Eqs. (|ll)6f) . I|1U7|) . (flXUp . 
pilfl . but nonzero values of Cs +Os are preferred. More 
precisely, in Solution I C,j g + ( g cannot be zero, while it 
is still consistent with zero at ~ 1.5a in Solution II (the 
poorly constrained orthogonal combination C,j g — Q g is 
consistent with zero in both cases). 
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